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Introducción 


Sin duda, el mayor acontecimiento astronómico de la década de los ochenta es 
la nueva visita del cometa Halley a la Tierra. Esta es su cita número 30 a la que 
puntualmente acude desde el año 239 a.C. 

El Halley se encontrará con dos grandes novedades científicas. Por primera 
vez el hombre saldrá al encuentro del cometa y hasta intentará entrar en contacto 
físico, íntimo, con él. El proyecto “Giotto” es la empresa internacional más apasio- 
nante en este aspecto, ya que tratará de que una nave espacial se acerque tanto 
al núcleo del cometa que casi todos los indicios apuntan a una destrucción 
irremediable (podrás observar este hecho en uno de los programas). Otras naves 
espaciales, rusas, japonesas y americanas, observarán al cometa desde el espacio. 

El otro acontecimiento no es menos importante. El cometa se acerca a un 
mundo nuevo, totalmente diferente al que vio en sus anteriores periplos. La gene- 
ración que habita actualmente el planeta Tierra le recibe amistosamente, sin recelos 
ni supercherías. Es una generación joven que no tiene miedo a acudir a su en- 
cuentro y que es capaz de observarle y estudiarle como un fenómeno de la natura- 
leza al que indudablemente respeta y admira y, hasta podríamos decir, tiene cariño. 

Esta nueva generación añade al conocimiento de la naturaleza astronómica del 
cometa la posibilidad de poder estudiarlo de forma extraordinariamente sofisticada 
desde su propia casa, con un instrumento cuya extensión masiva era inimaginable 
hace tan siquiera unos años: es el ordenador doméstico. 

La presente publicación surgió entre la efervescencia de la visita del Halley 
que ha afectado a los sectores más activos de los amantes de la astronomía y 
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de la ciencia. La idea nació y se desarrolló entre conferencias abiertas, reuniones 
de miembros de los grupos astronómicos, excursiones nocturnas para la observa- 
ción de “nuestro amigo Halley”, apariciones en los medios de difusión... En medio 
de ello destacaba la presencia de este nuevo elemento que facilitaba enormemente 
la comprensión de cada uno de los fenómenos explicados por los astrónomos: 
el microordenador. Para muchos dejó de ser un “juguete” para convertirse en un 
medio idóneo para la comprensión de las leyes que rigen el movimiento de los 
cuerpos celestes. 

Dos partes componen esta obra. La primera es una magnífica introducción al co- 
nocimiento de los fenómenos astronómicos destacando el papel de los cometas. Ha 
sido elaborada por Agustín Sánchez López, secretario de la Agrupación Astronó- 
mica Vizcaína, en un esfuerzo integrador que incluye los dos grandes avances: 
la intervención de la informática en el mundo de la astronomía y el apasionante 
relato del programa internacional para la recepción espacial del Halley. 

La segunda parte incluye los programas de ordenador por medio de los cuales 
puedes comprender diversos aspectos de las leyes que rigen el sistema solar 
y los cometas. Ha sido elaborada por miembros de la Agrupación Astronómica 
Salmantina. Los programas se presentan de forma que la comprensión de los con- 
ceptos es paulatina desde lo más simple hasta lo más complejo. Se ha cuidado con 
esmero la correlación entre los resultados en pantalla y la realidad astronómica. 


Salamanca, a 8 de diciembre de 1985 
KENNETH J. THOMSON 
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El Sol 


uestro Sistema 
Planetario 


Nuestro Sistema Planetario está constituido por los siguientes cuerpos: el Sol, 
nueve planetas, al menos 44 satélites y una infinidad de cuerpos mucho más pe- 
queños, que son: los asteroides, los meteoritos, los micrometeoritos y los cometas. 


El Sol es el cuerpo central del Sistema. En realidad, se trata de una estrella, 
más bien pequeña, dentro de la familia estelar, a pesar de lo cual representa el 
99,9 por 100 de la masa total del Sistema Planetario. Por ello, él es el centro de 
atracción gravitatoria del Sistema y el que regula el movimiento de los demás 
cuerpos, obligándoles a describir Órbitas a su alrededor y arrastrando a todos 
ellos en su movimiento por el espacio. 

El Sol es una inmensa esfera incandescente de, aproximadamente, 1.400.000 Km 
de diámetro, es decir, unas 109 veces mayor que el diámetro terrestre, y que 
posee una masa 330.000 veces mayor que la de la Tierra, encontrándose de 
nosotros a una distancia media de unos 150 millones de Km. 

También sabemos que funciona como un gigantesco reactor termonuclear, ob- 
teniendo su energía por la fusión del hidrógeno en helio, quemando en este pro- 
ceso la fantástica cifra de 600 millones de toneladas de hidrógeno por segundo. 

A pesar de tan impresionante consumo, lleva funcionando desde hace unos 
4.500 millones de años, y todavía le queda suficiente combustible nuclear para 
otros 4.500 millones de años más. 
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Los planetas 


En la actualidad se conocen nueve, de dimensiones y constitución variable. 
Por orden de proximidad al Sol son: Mercurio, Venus, Tierra, Marte, Júpiter, 
Saturno, Urano, Neptuno y Plutón. 

Atendiendo a la situación que ocupan en el espacio, en relación con la posición 
de la Tierra, los podemos clasificar en dos clases: planetas interiores y planetas 
exteriores. 

Los interiores, también llamados inferiores, son los que se encuentran en el 

interior de la órbita de la Tierra. Son solamente dos: Mercurio y Venus. Los ex- 
teriores, todos los demás, y se hallan fuera de la órbita de la Tierra. 

Si tenemos en cuenta algunas de sus características, como tamaño, densidad, 
tipos de atmósferas que poseen, velocidades de rotación, etc., los planetas admiten 
otra clasificación en dos tipos: los planetas telúricos o terrestres y los planetas 
jovianos o gigantes. 

En el primer grupo se encuentran Mercurio, Venus, Tierra y Marte. En el se- 
gundo, Júpiter, Saturno, Urano y Neptuno. En esta clasificación no entra Plutón, 
por no poseer las características comunes a los cuatro planetas jovianos, quizá 
por tratarse, como algunos investigadores piensan, de un antiguo satélite de Nep- 
tuno que fue arrancado de su órbita por la influencia gravitatoria del satélite 
Tritón. 


Los satélites 


Algunos de los planetas poseen, a su vez, otros cuerpos de menores dimen- 
siones que ellos orbitando a su alrededor. Son los que llamamos satélites. 

Por lo que conocemos hasta la fecha, Mercurio y Venus no poseen ninguno; 
la Tierra tiene uno, la Luna; Marte tiene dos, Júpiter posee dieciséis, Saturno 
diecisiete, Urano cinco, Neptuno dos y Plutón uno. En total, 44 satélites repar- 
tidos de forma irregular entre los nueve planetas. 


Los asteroides 


Además de los planetas y satélites, nuestro Sistema Planetario cuenta con unos 
50.000 cuerpos de dimensiones comprendidas entre 1 y 800 Km de diámetro. 
Son los que llamamos asteroides, la mayoría de los cuales se encuentran orbitando 
alrededor del Sol entre las órbitas de Marte y Júpiter, formando lo que llamamos 
el “cinturón de asteroides”. 


Los meteoritos y micrometeoritos 
El espacio interplanetario está inundado de cuerpos de reducidas dimensiones. 


Son los que llamamos meteoritos y micrometeoritos, de dimensiones compren- 
didas entre algunas micras y 1 Km. 
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En numerosas ocasiones, la Tierra, en su movimiento a través del espacio, se 
cruza con algunos de estos meteoritos, sometiéndoles a su influencia gravitatoria. 
Progresivamente los acelera, por la acción continua de esa fuerza, hasta obligarles 
a alcanzar velocidades que superan los 100.000 Km por hora. A estas velocidades, 
esos cuerpos, que a veces no son mayores que un grano de arena, llevan una 
enorme energía cinética acumulada, la cual se transforma en calor por rozamiento 
al penetrar en nuestra atmósfera, provocando la ionización de las moléculas del 
aire y volatilizando las partículas constituyentes del meteorito. 

Este efecto se manifiesta visualmente por esos fenómenos luminosos tan espec- 
taculares que llamamos meteoros y estrellas fugaces. 

A veces, algunos de estos cuerpos son tan grandes que no llegan a volatilizarse 
totalmente en nuestra atmósfera, haciendo entonces impacto en la superficie terres- 
tre, formando cráteres más o menos grandes, en función del tamaño del meteorito 
que cae. Entre los cráteres más espectaculares formados por un meteorito se en- 
cuentra el del Cañón del Diablo de Arizona (Estados Unidos), que tiene un diá- 
metro superior a 1 Km y una profundidad de unos 200 m. Se trata del resultado 
del impacto de un meteorito ocurrido hace unos 25.000 años. 

Los micrometeoritos son material pulverulento, en su mayoría de tamaño mi- 
croscópico, que bombardean sin cesar la atmósfera terrestre, asentándose lenta- 
mente en su superficie. 

Los micrometeoritos son tan numerosos que se estima se depositan sobre la 
superficie de la Tierra en cantidades del orden de las 100.000 toneladas al año. 


Los cometas 


Finalmente tenemos los cometas. Estos cuerpos celestes presentan un aspecto 
muy peculiar, distinguiéndose tres partes en su estructura morfológica, perfecta- 
mente diferenciadas: la zona más brillante, llamada “núcleo fotométrico” (para 
distinguirlo del verdadero núcleo, tan pequeño que resulta imposible observar 
desde la Tierra ni con el más potente telescopio); la nebulosidad que lo rodea, 
llamada “coma” o “cabellera”, y ese curioso apéndice llamado “cola”. 

Precisamente por presentar estas extrañas formas y por aparecer en el cielo 
siempre de una manera súbita e inesperada, los cometas han atemorizado al 
hombre desde la más remota antigiiedad, considerándolos heraldos de cataclismos 
y portadores de presagios funestos. 

En la actualidad, aunque los cometas siguen siendo los grandes desconocidos 
del Sistema Planetario, existen varias hipótesis fundadas en sólidos datos obtenidos 
por numerosas observaciones, que nos proporcionan gran información acerca de 
su naturaleza, su origen y del porqué de sus extrañas formas y sus movimientos. 


Distancias planetarias 
Con el fin de hacernos una correcta idea de la escala de nuestro Sistema 


Planetario, veamos cuáles son las distancias más destacables entre los diferentes 
miembros del Sistema. Por ejemplo, la distancia media del Sol a la Tierra es de 
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1. Núcleo verdadero.—2. Cola.—3. Núcleo fotométrico.—4. Cabellera o coma 


unos 150 millones de Km y la distancia media del Sol al planeta más alejado, 
Plutón, es de unos 6.000 millones de Km. 

Utilizando la “unidad astronómica de distancia” (u.a.), muy empleada en Astro-- 
nomía al referirse a las distancias interplanetarias, teniendo en cuenta que 1 u.a. 
es precisamente, la distancia del Sol a la Tierra, tendremos que el diámetro de la 
hipotética esfera sería de unas 80 u.a. 

Si empleamos el tiempo luz, es decir, la distancia recorrida por la ls en un 
determinado tiempo, sabiendo que la luz recorre, aproximadamente, 300.000 Km 
en un segundo, tendríamos un valor de unos 8 minutos luz para la distancia 
media de la Tierra al Sol y unas 6 horas luz para la distancia media del Sol 
a Plutón. 


Nuestra galaxia 


Nuestro Sol no es más que una de los 100.000 millones de estrellas que se 
agrupan entre sí para formar nuestra galaxia, llamada «Vía Láctea». 

La Vía Láctea es una galaxia típica, que tiene forma de disco con brazos 
espirales. Su diámetro es de unos 100.000 años luz y el espesor de su núcleo central 
de unos 20.000 años luz. 

El Sol, con su cortejo de planetas, satélites y cuerpos menores, ocupa el 
extremo de unos de esos brazos espirales, el llamado “brazo local o de Orión”, 
que se encuentra situado entre el “brazo de Perseo” y el “brazo de Sagitario”. 

Su distancia al núcleo es de unos 30.000 años luz, es decir, las tres quintas 
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partes del radio de la galaxia, y gira a su alrededor con una velocidad de casi 
1 millón de Km por hora, empleando unos 200 millones de años en dar una vuelta. 

Nuestra Vía Láctea se encuentra, a su vez, asociada con una veintena más de 
galaxias, formando el llamado “grupo local de galaxias”, que está formado por 
tres galaxias espirales, incluida la nuestra, varias elípticas y algunas irregulares. 

Entre las irregulares destacan las dos Nubes de Magallanes, visibles desde el 
hemisferio austral; se trata de dos galaxias satélites a la nuestra, situadas a unos 
200.000 años luz de distancia. 

Entre las espirales, la más importante es M-31, mejor conocida como galaxia 
de Andrómeda, que es la única visible a simple vista desde nuestro hemisferio. 
Esta galaxia se encuentra de nosotros a una distancia de unos 2 millones de años 
luz y posee un número de estrellas y una forma muy similar a la Vía Láctea. 


Distancias estelares y galácticas 


Hemos visto que el planeta más alejado del Sol, Plutón, se encuentra a 6 horas 
luz de aquél. Pues bien, para llegar a la estrella más cercana a nosotros, Pró- 
xima de Centauro, tenemos que dar un salto de 4,3 años luz, lo que nos da una 
idea del inmenso vacío existente en el espacio interestelar. 

Si queremos llegar al centro de la Vía Láctea, tendremos que dar un nuevo 
salto hasta los 30.000 años luz. Para alcanzar las galaxias más próximas necesita- 
ríamos 200.000 años luz y 2 millones para alcanzar la galaxia de Andrómeda, 
teniendo que aumentar estas distancias hasta valores de varios miles de millones 
de años luz para alcanzar los límites del universo actualmente conocido. 
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la naturaleza 
de los cometas 


A pesar de los grandes avances experimentados por la Astronomía en las 
últimas décadas, los cometas siguen siendo en la actualidad los grandes descono- 
cidos del Sistema Planetario. 

Todo lo que la Astrofísica moderna conoce sobre los cometas debe dividirse 
en dos categorías: lo que podemos considerar como certeza y lo que son simples 
hipótesis, que todavía tienen que ser confirmadas por medios directos. 

Entre los conocimientos ciertos se encuentran todos los aspectos relacionados 
con su dinámica y su constitución externa. 

En la categoría de las hipótesis están su constitución interna, su proceden- 
cia y su origen. 


Aspecto de los cometas 


Según dijimos antes, el aspecto que presentan los cometas es el de un punto 
más o menos brillante, el “núcleo fotométrico”, que resulta perfectamente visible 
a través del telescopio, rodeado de una débil nebulosidad llamada “coma” o “ca- 
bellera”, acompañada de un curioso apéndice llamado «cola», que siempre se 
encuentra orientada en dirección contraria al Sol, por los motivos que más ade- 
lante veremos. 

El tamaño de las cabelleras es muy variable, encontrándose la mayoría de 
ellas entre los 15.000 y 250.000 Km de diámetro. 
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La mayor cabellera observada fue la del cometa de 1811, que alcanzó un diá- 
metro de 1.800.000 Km, es decir, mayor que el propio Sol. 


Mucho más asombrosa aún es la extensión de sus colas, que suelen alcanzar 
varios millones de kilómetros de longitud. 

La cola del cometa Halley, en su presentación de 1910, fue de 110 millones 
de Km. La cola más larga observada fue la del cometa de 1843, que llegó a alcan- 


zar un valor de 320 millones de Km. Es decir, muy superior a la distancia que 
separa el planeta Marte del Sol. 


Cola de plasma 
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Movimiento del cometa 
Coma 


Núcleo fotométrico 


Hacia el Sol 


Aspecto global de un cometa 


¿Qué son los cometas? 


Para contestar a esta fundamental pregunta, hemos de recurrir a la hipótesis 
formulada en el año 1950 por el astrónomo norteamericano Fred Whipple. 

Para Whipple, los cometas simplemente son cuerpos muy pequeños, de sólo 
unos pocos kilómetros de diámetro (entre 5 y 10 Km) y con una forma más o 
menos irregular, pero tendiendo a la esférica, que están formados por un conglo- 
merado de rocas, polvo y hielos. A este modelo del núcleo de un cometa se le 
conoce como el modelo de la “bola de nieve sucia”. 


¿Cómo estas pequeñas bolas se presentan 
de forma tan espectacular? 


¿ 


La explicación es muy sencilla tomando el modelo de núcleo cometario de 
Whipple. Cuando esa bola de nieve sucia, que flota a través del espacio, se apro- 
xima al Sol, sus capas externas constituidas por partículas de polvo y hielos de 
varios compuestos químicos como el agua, el amoníaco, el metano, el anhídrido 
carbónico, el ácido cianhídrico, etc., son calentados por la radiación solar, lo que 
provoca la sublimación de los hielos que se convierten en gas. 

Estos/ gases sublimados se expanden en el espacio interplanetario, arrastrando 
consigo a las partículas de polvo, que luchan contra la débil gravedad del pequeño 
núcleo, para formar una nube gaseosa y polvorienta que rodea al núcleo. Esta 
nube de gas y polvo constituye la cabellera del cometa, la cual brilla, funda- 
mentalmente, al difundir la luz que reciben del Sol. 


Cómo se forma la gigantesca cola? 


A finales del siglo xIX, el matemático y físico escocés James Clerk Maxwell 
demostró teóricamente que la luz ejerce una presión sobre los cuerpos que en- 
cuentra en su camino. Ello es debido a la naturaleza corpuscular de la luz, co- 
nociéndose este efecto como “presión de la radiación”. 

Esta presión de la radiación solar arrastra al polvo de la cabellera, formando 
las características colas cometarias, las cuales, por su propia naturaleza, necesaria- 
mente tienen que estar orientadas en dirección contraria al Sol. 

Como la presión de la radiación solar decrece con el cuadrado de la distancia, 
se comprende fácilmente que a medida que el cometa se va acercando al Sol, 
las colas serán cada vez mayores, alcanzando su máximo esplendor en el momento 
del paso del cometa por el perihelio. 

El primero en sugerir este mecanismo de formación de las colas cometarias fue 
el físico irlandés George Francis Fitzgerald, en el año 1882, pero erróneamente 
consideró que en ese mecanismo estaban incluidas las moléculas de los gases, 
dándose cuenta posteriormente que la presión de la radiación no ejerce una 
fuerza suficiente sobre las moléculas que permitan explicar la formación de las 
colas. 

Fue el químico sueco Arrhenius quien, en el año 1900, emitió la hipótesis de 
que la materia repelida por la presión de la radiación no son moléculas de 
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Orientaciones y dimensión de la cola de un cometa en su paso alrededor de una estrella 


gas, sino que se trata solamente de partículas pequeñísimas de polvo. Las molé- 
culas, al ser más pequeñas que la longitud de onda de la luz, en general, no 
interaccionan con los fotones de la luz y, por tanto, no pueden ser rechazadas 
por la radiación solar. 

Por el mecanismo descrito, queda perfectamente explicada la formación de las 
colas de polvo de los cometas. Estas colas se conocen como “colas del Tipo IT”, 
para distinguirlas de otras muy diferentes que también se observan en la mayoría 
de los cometas, llamadas “colas de plasma” o “colas del Tipo TI”. 

Las colas de polvo son curvadas, anchas y relativamente cortas, pues en general 
no sobrepasan los 10 millones de Km de longitud. Presentan una tonalidad ama- 
rillenta-rojiza, característica de la luz solar difundida por las partículas de polvo 
que la forman. Estas partículas de polvo tienen unas dimensiones comprendidas 
entre 0,2 y 5 micras. 

Por el contrario, las colas del Tipo 1, o colas de plasma, son rectilíneas, 
mucho más largas que las de polvo y presentan una tonalidad azulada. 

Estas colas de plasma están constituidas por gas ionizado, debiendo su color 
azulado a la emisión por fluorescencia del ion CO*, que emite con una longitud de 
onda de 4.200 amstrong. 

Pero, ¿cuál es el mecanismo de formación de estas colas de plasma? Esta 
pregunta intrigó durante mucho tiempo a los astrónomos y no pudo ser correc- 
tamente contestada hasta el año 1951. 

Fue el astrofísico alemán Ludwig Biermann quien sugirió por primera vez en 
1951 que las colas de plasma se debían a la acción del «viento solar» sobre las 
moléculas de la cabellera. Este viento solar está formado por un flujo de partículas 
subatómicas, esencialmente protones, electrones y núcleos de helio, junto con otras 


pequeñas cantidades de iones de elementos más pesados, como el oxígeno y el 
carbono. 

El Sol está expulsando continuamente un millón de toneladas de gas por se- 
gundo, a una velocidad de unos 400 Km por segundo. Las partículas de alta 
energía de este viento solar junto con la radiación ultravioleta de la luz, producen 
la ionización de las moléculas de gas de la cabellera, formándose iones CO*, 
CO y N3 y radicales OH* y CH*, los cuales son arrastrados por los campos 
magnéticos asociados al viento solar, acelerando dichos iones hasta velocidades 
de varias decenas de kilómetros por segundo. 


¿De dónde proceden los cometas? 


Nuevamente nos encontramos con otro interrogante para cuya respuesta hemos 
de recurrir a otra hipótesis. En esta ocasión fue el astrónomo holandés Jan 
Hendrik Oort quien, en el año 1950, formuló su hipótesis de que deben existir del or- 
den de 200.000 millones de núcleos cometarios, de “bolas de nieve sucia” según el 
modelo de Whipple, hibernando en la profunda gelidez del espacio interestelar, 
en los confines del Sistema Planetario, es decir, en los límites de la esfera de in- 
fluencia gravitatoria del Sol. 

Estos límites se encuentran, aproximadamente, a un año luz de distancia del 
Sol, donde los núcleos cometarios se encuentran formando una especie de nube, 
conocida como la “nube de Oort”,en honor de su descubridor. 

Todos estos núcleos se encuentran sometidos a la débil influencia gravitatoria 
que el Sol ejerce a tan enorme distancia, describiendo órbitas elípticas de muy 
baja excentricidad, es decir, casi circulares. 

A pesar de existir tan gran cantidad de núcleos cometarios en la Nube de 
Oort, al ser estos núcleos cuerpos tan pequeños, la masa total de la Nube es 
sensiblemente inferior a la masa de la Tierra, estimándose un valor de sólo la dé- 
cima parte de ésta. 

Unos cuerpos tan poco masivos y situados a tan lejana distancia del Sol, nece- 
sariamente tienen que resultar muy afectados por la influencia gravitatoria de las 
estrellas vecinas. 

Las perturbaciones gravitacionales que introducen dichas estrellas en sus movi- 
mientos relativos respecto al Sol, son las causantes de que en ocasiones sean 
arrancados de la Nube un número determinado de núcleos cometarios, algunos 
de los cuales son lanzados al espacio exterior, otros quedan atrapados alrededor 
de la estrella y otros, en cambio, son desviados en sus trayectorias para adentrarse 
en el interior del Sistema Planetario. 

A medida que se van aproximando al Sol irán formando su cabellera y sus 
colas por el mecanismo que antes hemos descrito, siendo modificadas sus órbitas, 
que inicialmente eran parabólicas, por otras elípticas, debido a la influencia gravi- 
tatoria de los planetas, convirtiéndose así en cometas periódicos. 

Para Oort, el origen de esta nube de cometas, hemos de buscarlo en el propio 
origen del Sistema Solar, a partir de la nube de gas y polvo primitiva, siendo los 
núcleos cometarios, por tanto, los restos fosilizados de esa nube. 
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Los cometas 
en la historia 
de la humanidad 


Hemos visto en el capítulo anterior cuáles son los conocimientos que la ciencia 
actual posee sobre los cometas. 

Pero, ¿cuáles han sido los pasos más destacables que la humanidad ha dado 
en el curso de la historia para llegar al conocimiento actual? 

Para los filósofos de la escuela pitagórica, en el siglo vI a.C., los cometas 
eran considerados como planetas que se elevaban muy poco por encima del hori- 
zonte, por lo cual pasaban inadvertidos la mayoría de las veces. 

Para Anaxágoras y Demócrito eran el resultado de choques de cuerpos celestes, 
y para Diógenes y Apolonio se trataban de estrellas. 


El pensamiento de Aristóteles 


El pensamiento de Aristóteles en relación con los cometas tiene una importancia 
trascendental en la historia, ya que sus ideas, a pesar de ser equivocadas, preva- 
lecieron hasta el siglo XVII, al igual que ocurriría con otras muchas cuestiones re- 
lacionadas con la ciencia en general y con la Astronomía en particular, cegando 
así a la civilización occidental durante más de 2.000 años. 

Para Aristóteles, el universo estaba dividido en dos regiones distintas: la región 
«sublunar», situada debajo de la Luna, donde todo era cambio y decadencia, y 
la región más alejada o «supralunar», donde todo era inmutable y eterno. 
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En sus Meteorológicas, Aristóteles escribiría que los cometas eran llamaradas 
de emanaciones singulares de nuestra atmósfera, situados en la región sublunar, que 
se encendían de vez en cuando y ardían durante semanas o meses, para luego, 
poco a poco, irse extinguiendo. 


La genial intuición de Séneca 


Séneca, filósofo cordobés del siglo 1, expuso ideas más acertadas. Para él, 
los cometas debían ser incluidos en la categoría de los planetas y, por tanto, 
como éstos, deberían tener un curso regular perfectamente calculable. 

Así pues, escribiría: «Sería indispensable poseer el catálogo de todos los cometas 
que aparecieron en el pasado. Todavía no podemos, debido a su rareza, conocer 
su camino, ni saber si su regreso es periódico y si un orden establecido los trae 
de nuevo en una fecha concreta.» 

La genial intuición del filósofo cordobés, al aludir a la posible periodicidad de 
los cometas, es verdaderamente sorprendente, pues es necesario esperar diecisiete 
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El cometa Halley en 1910, 


siglos para verla confirmada, tras la predicción que efectuaría Edmund Halley, 
según veremos más adelante. 

Séneca profetizó lo siguiente: «Día vendrá en que lo que hoy es misterioso 
será puesto a la luz. Nacerá un hombre que hará ver en qué parte del cielo van 
los cometas, por qué se diferencian tanto de los demás planetas y cuál es su mag- 
nitud y su naturaleza.» 


La primera medición de la distancia de un cometa 


Pero el primer paso histórico intentando demostrar el error de la concepción 
aristotélica, lo daría el astrónomo alemán Johann Muller, mejor conocido por el 
seudónimo de Regiomontano. 

Regiomontano estaba convencido de que los cometas eran astros y no meteoros, 
como sostenían los partidarios de Aristóteles. Para demostrar su idea, intentó me- 
dir la distancia a la que se encontraba un cometa muy luminoso, que apareció 
en el año 1456 (como veremos, se trataba del cometa Halley). 
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Para ello, utilizó el método del paralaje, que consiste en medir el desplaza- 
miento angular que sufre un objeto, respecto al fondo de las estrellas fijas, reali- 
zando la medición desde dos puntos separados entre sí una distancia conocida. 
Las mediciones que efectuó no tuvieron la suficiente precisión y obtuvo un valor 
de sólo nueve veces el radio de la Tierra, es decir, aproximadamente la séptima 
parte de la distancia de la Tierra a la Luna. 

Con los datos obtenidos, no podía refutarse la hipótesis de Aristóteles, pues 
el cometa parecía encontrarse en la región sublunar. 


Las colas planeatarias sometidas a las leyes 
de la naturaleza 


Setenta y cinco años después, dos astrónomos, el italiano Jerónimo Fracastoro y 
el alemán Pedro Apiano, trabajando por separado, observaron al cometa Halley en 
su presentación del año 1531, llegando ambos a la misma conclusión: la cola del 
cometa apuntaba siempre en dirección contraria al Sol. Por primera vez se 
demostraba que los cometas seguían leyes fijas. 


Copérnico: una revolución en el orden establecido 


El astrónomo polaco Nicolás Copérnico, en el año 1543, desencadenó la mayor 
revolución científica y cultural de la historia de la humanidad, al afirmar que era el 
Sol y no la Tierra quien ocupaba el centro del Sistema Planetario. 

Copérnico sentaba de esta manera las bases de un sistema muchísimo más 
sencillo que el desarrollado por el astrónomo alejandrino Ptolomeo, que era el gene- 
ralmente aceptado en esa época. Pero la nueva concepción de Copérnico trastocaba 
todo el orden establecido y, además, parecía estar en contradicción con algunos 
pasajes de las Sagradas Escrituras. Todo ello, hacía que esta afirmación fuera con- 
siderada herética por los encargados de preservar las ideas antiguas. 

Por defender esta opinión, entre otras que hoy consideramos naturales, Gior- 
dano Bruno fue quemado en la hoguera por la Inquisición en el año 1600 y 
Galileo Galilei sufriría persecución, encarcelamiento y tortura, siendo obligado 
finalmente a retractarse en el año 1633, 


Tycho Brahe: un excelente observador 
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Avanzando en la historia, nos encontramos con otro gran personaje: el danés 
Tycho Brahe, famoso por ser el astrónomo de su época que efectuaba con mayor 
precisión las medidas astronómicas. 

En el año 1572 ocurrió un fenómeno en el firmamento decisivo para asestar 
un golpe definitivo a la idea aristotélica del universo. En este año, aparecía en el 
cielo una estrella nueva, tan luminosa como el propio Júpiter, en un punto del cielo 
donde antes no se había observado nada. 


Kepler: 


Tycho Brahe no quería creer lo que veía con sus propios ojos. Para los sabios 
europeos de la época, tenía que tratarse de un fenómeno de la región sublunar, 
la supuesta «región de los cambios». 

Tras laboriosas mediciones, las conclusiones finalmente llegaron: sin lugar a du- 
das, esa estrella se encontraba muchísimo más alejada que la Luna, incluso, mucho 
más alejada que todos los planetas conocidos en la época. Aquella «nueva» estrella 
se encontraba situada en la octava esfera de Ptolomeo, es decir, entre el resto de 
las estrellas fijas, en la región considerada como «inmutable». 

Hoy sabemos que lo que Tycho vio no fue el nacimiento de una nueva estrella, 
sino el último estertor de una estrella muy masiva de nuestra galaxia, que había 
llegado al fin de su vida, brillando esplendorosamente en sus últimos momentos, 
cual si de un Ave Fénix se tratara. 

El 13 de noviembre de 1577, cuando Tycho se dirigía a un estanque para pes- 
car, descubrió un cometa extraordinariamente brillante, tanto, que incluso se veía 
a la luz del Sol. 

Dispuesto a resolver definitivamente el problema, utilizó el mismo método 
del paralaje empleado anteriormente por Regiomontano, pero con una precisión 
jamás conseguida anteriormente. Los datos obtenidos eran definitivos. El cometa 
en cuestión se encontraba, al menos, seis veces más alejado que la Luna. 

Esta era la prueba definitiva contra las teorías de Aristóteles. Los cometas no 
eran fenómenos atmosféricos, sino cuerpos celestes de pleno derecho. 


las leyes del movimiento planetario 


Al astrónomo alemán Johann Kepler, colaborador de Tycho, le correspondía 
inaugurar, en el año 1609, la moderna astronomía planetaria al descubrir sus in- 
mortales leyes del movimiento de los planetas. 

La primera de estas leyes indica que los planetas describen órbitas elípticas al- 
rededor del Sol, que ocupa uno de los focos. 

La segunda explica que los planetas se mueven más rápidamente en su órbita 
cuanto más cerca se encuentre del Sol, con una relación perfectamente determinada. 

La tercera relaciona los cuadrados de los tiempos empleados en describir las 
órbitas, con los cubos de los semiejes mayores de ellas. 

Con estas leyes, tenemos la explicación cinemática de los extraños movimientos 
que realizan los planetas contra el fondo de las estrellas. 

Pero, y los cometas, ¿cómo se mueven? Kepler observó un cometa que apa- 
reció en el año 1607. Se trataba nuevamente del cometa Halley. Posteriormente 
observó otro distinto en el año 1618. Según sus observaciones, estos cuerpos no 
seguían las leyes por él enunciadas. Para Kepler, los cometas viajan eternamente 
a una velocidad regular y en línea recta, cruzando el universo infinito. 


Isaac Newton 


Isaac Newton fue muy probablemente el genio más grande de la historia de la 
humanidad. El sería el hombre que profetizara Séneca dieciséis siglos antes. 
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¿Cuáles fueron las aportaciones de Newton a la ciencia? Sin lugar a dudas, 
la más importante de todas ellas fue el descubrimiento de su famosa ley de la 
Gravitación Universal, publicada en el año 1687, en su libro Principios matemáti- 
cos de la filosofía natural. 

Según esta ley, todos los cuerpos del universo se atraen con una fuerza que es 
directamente proporcional a sus masas e inversamente proporcional al cuadrado 
de sus distancias. 

Esta es precisamente la fuerza que provoca que los planetas describan órbitas 
elípticas alrededor del Sol, según había descubierto Kepler setenta y ocho años 
antes. 

Ya tenemos la explicación dinámica que nos faltaba: Los planetas son, en reali- 
dad, masas en movimiento bajo la influencia del campo gravitatorio del Sol. 

Y los cometas, ¿estaban sometidos a la misma fuerza? Para Newton la contes- 
tación era evidente. No existía ninguna razón para que los cometas se escaparan 
a la influencia de esta fuerza. 

Entonces, ¿se movían siguiendo órbitas elípticas alrededor del Sol como los pla- 
netas? ¿No seguían un movimiento rectilíneo, como argumentaba Kepler? Newton 
observó un cometa en el año 1680 y su movimiento parecía corresponderse con 
una parábola que tuviera como foco el Sol. Pero él estaba convencido de que las 
trayectorias de los cometas tenían que ser elipses, como las de los demás planetas. 
Eso sí, elipses muy alargadas, las cuales serían muy difíciles de distinguir de las 
parábolas. 


Edmund Halley 
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Nos encontramos ya con nuestro protagonista principal, el astrónomo inglés 
Edmund Halley, contemporáneo y amigo de Newton. Halley estaba empeñado en 
situar a los cometas en su verdadero lugar dentro del orden cósmico y, por con- 
siguiente, hacer desaparecer todas las especulaciones y temores que estos cuerpos 
provocaban. 

También él observó el cometa de 1680. Al igual que Newton, comprobó que 
se había estado acercando al Sol en noviembre de ese año, para después perderse 
por encima del Sol en diciembre, siguiendo una dirección casi completamente 
opuesta. Esto era una prueba evidente de que tenía que haber cambiado su curso. 

Según Newton explicó a Halley, la trayectoria tenía que ser aproximadamente 
una parábola. 

Le correspondía a Halley perfeccionar la obra de Newton demostrando que los 
cometas tenían una órbita elíptica y que, por tanto, son astros periódicos. 

Halley observó un nuevo cometa en el año 1682. Este sería el cometa que le 
daría la inmortalidad. 

Dispuesto a resolver el enigma de los movimientos de estos cuerpos, abordó la 
ardua tarea de calcular las órbitas de veinticuatro cometas de los que poseía un 
número suficiente de datos. Para ello utilizó el método matemático elaborado por 
el propio Newton. 

- Las conclusiones no se hicieron esperar: En 1705 publica su memoria Astro- 


nomica Cometicae Synopsis, y escribe lo siguiente: «Muchas cosas me inducen a 
pensar que el cometa del año 1531 observado por Apiano es el mismo que fue des- 
cubierto por Kepler en 1607 y que personalmente he observado yo en su regreso 
de 1682.» 

Para Halley, por tanto, los tres cometas eran el mismo, que describía una 
órbita elíptica muy excéntrica y que se acercaba a la Tierra cada 75 6 76 años. 
Estos tres cometas se movían en sentido retrógrado (contrario al de todos los 
planetas del Sistema Solar), al igual que lo hiciera el cometa del año 1456; 
por tanto, tenía que tratarse del mismo cometa. 

Finalmente, Halley se atrevió a realizar su famosa predicción: «Puedo confia- 
damente afirmar que este cometa regresará en el año 1758. Si esta predicción se 
cumple, no habrá ya razón alguna para dudar de que regresarán otros cometas.» 

Halley falleció en 1742, a la edad de ochenta y seis años, dieciséis años antes 
de la fecha prevista. Cumpliéndose su predicción en la Navidad del año 1758, un 
astrónomo aficionado alemán, Johann Georg Palitzsch, anunciaba que el tan 
esperado cometa se encontraba ya al alcance de los telescopios de la época. Desde 
entonces, al cometa se le llama Halley, convirtiéndose así en el más famoso de 
todos los cometas periódicos. El cometa volvería a hacer su aparición en 1835 y 


en 1910, estando ya próximo su siguiente paso por el perihelio, el 9 de febrero 
de 1986. 
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Sucesos más 
destacables del 
cometa Halley en 
su retorno de 1986 


Sucesos antes del perihelio 


— 5 de mayo de 1948.—El cometa estaba situado en el afelio, es decir, en el punto 
de su órbita más alejado del Sol, a 35,3 u.a. (5.250 millones de Km), o sea, 
más alejado que Neptuno. La velocidad del cometa en ese punto era de sólo 
0,9 Km por segundo (3.240 Km por hora). 


— 16 de octubre de 1982.—Cuando se encontraba más alejado que Saturno (a 
11,1 u.a.) fue redescubierto por los astrónomos norteamericanos Jewitt y Da- 
nielson, por medio del telescopio de 5 m de diámetro de Monte Palomar y una 
cámara fotográfica electrónica ultrasensible (cámara CCD). 


— 10 de enero de 1985.—El cometa cruzó la órbita de Júpiter. 


— 24 de octubre de 1985.—Máximo de la lluvia de meteoritos «Las Oriónidas», 
asociadas al cometa Halley por tratarse de restos de materia expulsados por el 
núcleo del cometa, que se encuentran situados en el espacio interplanetario, 
siguiendo la misma órbita que el cometa. 


— 10 de noviembre de 1985.—Se encontrará en la máxima declinación Norte: 
+22%14". 


— 18 de noviembre de 1985.—En oposición con el Sol. Podrá ser observado durante 
toda la noche con unos simples prismáticos. Su ángulo de fase (ángulo Sol- 
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Cometa- Tierra) será de sólo 1,4”, el mínimo en esta presentación. La cola, por 
tanto, estará dirigida en la dirección de nuestra visual, resultando imposible 
su Observación. 


— 27 de noviembre de 1985.—Se encontrará a la mínima distancia de la Tierra 
antes del perihelio, a sólo 0,62 u.a. (unos 93 millones de Km). 


— 28 de noviembre de 1985.—Cruzará la órbita de Marte. 


— 25 de diciembre de 1985.—Máximo ángulo de fase antes del perihelio, aproxi- 
madamente 55”. Será un buen momento para estudiar la estructura de la cola. 
Probablemente comenzará a ser perceptible a simple vista de forma muy débil. 


— 1 de enero de 1986.—Cruzará la órbita terrestre. 
— 21 de enero de 1986.—Cruzará la órbita de Venus. 
— 6 de febrero de 1986.—Se encontrará en conjunción con el Sol. 


—-9 de febrero de 1986.—El cometa alcanzará el perihelio, el punto de máximo 
acercamiento al Sol. Se encontrará a 0,59 u.a. del Sol (88,3 millones de Km) y 
a 1,55 u.a. de la Tierra (231,9 millones de Km). Rebasará el perihelio a una 
velocidad de 54,5 Km por segundo (196.200 Km por hora). 


Sucesos después del perihelio 


— 1 de marzo de 1986.—Cruzará la órbita de Venus. Ahora en sentido descen- 
dente, contrario al que realizó antes del perihelio. 


— 8 de marzo de 1986.—La sonda japonesa Planet-A se acercará a unos 100.000 Km 
del cometa. 


— 9 de marzo de 1986.—La sonda japonesa MS-T5 cruzará la cola de plasma del 
cometa. 


— 10 de marzo de 1986.—La sonda soviética Vega 1 se acercará a 10.000 Km 
del cometa. 


— 13 de marzo de 1986.—La sonda europea Giotto se acercará a la increíble dis- 
tancia de 500 Km del núcleo del cometa. 


— 18 de marzo de 1986.—La sonda soviética Vega 2 se acercará a 7.000 Km del 
cometa. 


— 20 de marzo de 1986.—Máximo ángulo de fase en esta presentación. Alcanzará 
un valor de 66”, por lo que se tendrán muy buenas condiciones para el estudio 
de la estructura de la cola. 


— 21 de marzo de 1986.—El cometa cruzará la órbita terrestre, alcanzando por 
estas fechas su máximo brillo (magnitud aproximada 4). 
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—9 de abril de 1986.—Tendrá la máxima declinación Sur: —47%23'. Las condi- 
ciones de observación para los habitantes del Hemisferio Sur serán excelentes. 


— 11 de abril de 1986.—Minima distancia del cometa a la Tierra. Se encontrará 
a 0,42 u.a. (unos 63 millones de Km). 


— 17 de abril de 1986.—Se encontrará en oposición con el Sol, 
— 23 de abril de 1986.—Cruzará la órbita de Marte. 
— 25 de abril de 1986.—Dejará de ser observable a simple vista. 


— 2 de junio de 1986.—Dejará de ser visible a través de los pequeños telescopios 
de aficionado, pues comenzará a tener una magnitud más débil de 10. 
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Recomendaciones 
para la 
observación 

del cometa 


Con el fin de facilitar las observaciones del cometa, citamos a continuación las 
recomendaciones más importantes a tener en cuenta. 


Telescopios 


Emplear telescopios de corta distancia focal y máximo diámetro de abertura de 
objetivo posible, siempre con la mínima potencia que permita el telescopio, salvo 
cuando se desee realizar una observación en detalle del núcleo. 


Prismáticos 


Utilizar prismáticos cuyo valor de la pupila de salida esté comprendida entre 
7 y 8. Este dato se obtiene fácilmente, pues sólo tenemos que dividir el diámetro 
del objetivo en milímetros por los aumentos que proporcionan los prismáticos. 
Por ejemplo, unos prismáticos de 7 x 50 (muy económicos y apropiados para la 
observación de cometas) daría un valor de pupila de salida de 7,14. Todavía mejor, 
serían los de 11 x 80 (pupila de salida 7,27), pues al tener más diámetro de objetivo 
se pueden observar objetos mucho más débiles. Su problema es que son excesiva- 
mente caros. 
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Nunca utilizar prismáticos de elevados aumentos, pues además de necesitar 
imprescindiblemente un trípode para evitar el movimiento de la imagen, dan una 
pupila de salida menor, para el mismo diámetro de objetivo. Por ejemplo, un 
equipo de 16 x 50, tiene una pupila de salida de sólo 3,12, valor excesivamente 
bajo. 

Para elegir el lugar de emplazamiento y el momento más oportuno, deberá 
tenerse en cuenta lo siguiente: 


1.2 Huir de la iluminación y contaminación atmosférica de los núcleos urbanos. 
Hemos de tener en cuenta que la luz parásita de las ciudades, al difundirse 
en el polvo atmosférico, da lugar a un cielo con una cierta iluminación, 
aun en plena noche. Por ello, las observaciones deberemos realizarlas 
siempre en pleno campo. 


2.” Elegir el punto de observación de forma que, en la dirección en que se 
encuentre el cometa, no exista ninguna población, aunque ésta sea pequeña 
y aunque se halle muy alejada del observador, ya que su luz parásita se 
difunde en el horizonte impidiéndonos una correcta observación. 


3. Dado que en las fechas en que el cometa será más brillante, se va a 
encontrar muy bajo y, por tanto, cercano al horizonte, deberemos elegir 
un lugar de observación desde el que se domine un horizonte lejano y des- 
pejado de montes, pueblos, etc. 


4. Elegir noches en las que la Luna no se encuentre sobre el horizonte en el 
momento de la observación, pues su resplandor nos dificultaría la obser- 
vación del cometa. 


5.2 Finalmente, hemos de tener siempre presente que nuestros ojos necesitan 
un tiempo de adaptación a la oscuridad. Después de que transcurran 
algunos minutos en total oscuridad, nos daremos cuenta que somos capaces 
de distinguir los más finos detalles del cometa, detalles que en un principio 
nos resultaban imperceptibles. 


Las condiciones en que vamos a poder ver el cometa entre los meses de no- 
viembre de 1985 y mayo de 1986, son las siguientes: 


Mes de noviembre 
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Comenzará a ser visible a través de pequeños telescopios, presentándose como 
un punto brillante (el núcleo fotométrico) rodeado de una pequeña nebulosidad 
(la cabellera). La cola no podrá ser observada durante este mes, por encontrarse 
casi coincidente con la dirección de nuestra visual desde la Tierra (mínimo ángulo 
de fase de 1,4” el día 18). 

Durante todo el mes irá moviéndose a través de la constelación de Taurus, 
pasando muy próximo al cúmulo abierto de las Pléyades a mediados de mes, 
para terminar en Piscis a finales. 

A primeros de mes aparecerá por el horizonte Este al atardecer, poco después 
de la puesta del Sol. A mediados de mes aparecerá justo en el mismo momento de 


la puesta del Sol y un poco antes a finales. Por ello, se encontrará a una buena 
altura sobre el horizonte Sur a medianoche durante todo el mes. 

El 27 de noviembre será la fecha en que se aproximará más a la Tierra antes 
del perihelio. Se encontrará a unos 93 millones de Km. 


Mes de diciembre 


Comenzará a ser observable por medio de prismáticos desde primeros de mes 
y es probable que las personas con una visión muy buena puedan empezar a obser- 
varlo a simple vista (de forma muy débil) por Navidad. En estas fechas tendremos 
unas buenas condiciones de visibilidad transversal (ángulo de fase 557), por l lo que 
podrá observarse bastante bien la cola del cometa. 

Se moverá a través de la constelación de Piscis, terminando en Acuario a 
finales de mes. 

Poco después de la puesta del Sol, lo tendremos sobre el horizonte Sur 
durante todo el mes, aunque perdiendo cada vez más altura (63” para Madrid el 
día 1 y 48” el día 31). 


Mes de enero 


Seguirá siendo visible a simple vista, de forma muy débil, durante la primera 
quincena del mes, moviéndose a través de la constelación de Acuario. 


ESTE it 


1400 110 120 130 44h 150 460 470 480 490 200 2to 220 230 240 250 260 270 280 


"Posiciones del cometa Halley: (D En enero, 70-90 minutos tras la puesta del Sol, moviéndose por Piscis y 

Acuario y perdiendo altura sobre el horizonte. Y En marzo y abril, meses de mayor esplendor, 70-90 minutos 

antes de salir el Sol, en dirección Sur-Este, atravesando Sagitario (marzo). YY) En la segunda quincena de abril, 
70-90 minutos tras la puesta del Sol, ganando altura sobre el horizonte y alejándose de nosotros 
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Poco después de la puesta del Sol lo tendremos sobre el horizonte Sur-Oeste, 
aproximándose cada día más y más al Sol, hasta quedar sumergido en su resplandor 
en la segunda quincena del mes, dejando entonces de ser visible por impedírnoslo 
la luz del Sol. 


Mes de febrero 


Pasará por el perihelio el día 9, pero no tendremos ocasión de verlo durante 
prácticamente todo el mes, por encontrarse casi alineado con el Sol, quedando, 
por tanto, oculto en su resplandor. 


Meses de marzo y abril 


Serán los meses de mayor esplendor del cometa en esta presentación, pero los 
habitantes del Hemisferio Norte lo veremos muy próximo al horizonte, lo que 
restará espectacularidad al fenómeno. Quien tenga programado un viaje al He- 
misferio Sur, es el momento de realizarlo. La experiencia, sin duda, resultará 
inolvidable. Ciudades como Melburne, Santiago de Chile, Buenos Aires, Mon- 
tevideo o Johannesburgo resultarán óptimas para una excelente observación. 

Los que no podamos viajar a esos lugares, tendremos que conformarnos 
con una visión menos espectacular, para lo cual deberemos madrugar durante 
el mes de marzo para poder observarlo hacia el Sur-Este, unas horas antes de la 
salida del Sol. Lo encontraremos muy próximo al horizonte, claramente perceptible 
a simple vista, atravesando la constelación de Sagitario. 

El 20 de marzo se darán las mejores condiciones de visibilidad transversal 
(ángulo de fase 66”), por lo que la observación de la cola resultará óptima. 

El 11 de abril tendrá lugar la segunda aproximación a la Tierra, acercándose 
aún más que lo hiciera el 27 de noviembre. Lo tendremos a unos 63 millones 
de Km, pero no será observable desde prácticamente ningún punto de España 
por encontrarse casi en el horizonte. Esta circunstancia se dará durante toda la 
quincena de abril. 

Comenzaremos a verlo nuevamente en la segunda quincena de abril, pero ahora 
por la tarde, poco después de la puesta del Sol, hacia el Este, cada día un poco 
más elevado sobre el horizonte. 

Dejará de ser observable a simple vista a finales de abril. 


Mes de mayo 


Solamente podemos seguir observándolo con prismáticos durante la primera 
quincena y con pequeños telescopios durante la segunda, despidiéndonos del co- 
meta hasta el año 2062, en que nuestros hijos o nietos volverán a contemplarlo, 
quizá, en mejores condiciones que las que nos ha sido dado observarlo a nosotros. 
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La Informática 
y la Astronomía 


Desde que en la década de los cuarenta Von Newman consiguiera configurar la 
arquitectura básica de los ordenadores, la Astronomía, al igual que otras ramas de 
la ciencia, pronto empezaría a beneficiarse de las múltiples posibilidades que los 
nuevos equipos ponían a su alcance. 

Algunos de los enormes ordenadores que hoy pueden ser contemplados en el 
Museo del Ordenador de la ciudad de Boston, guardan una estrecha relación 
con las primeras aplicaciones de estos equipos en el campo de la ciencia. 

El famoso Whirlwind, desarrollado en 1950 por la compañía Digital Equip- 
ment Corporation, fue utilizado en el Instituto de Tecnología de Massachusetts 
(MIT) para fines científicos. El Whirlwind, fue el primero de los ordenadores al 
que se le equipó con memoria de núcleos magnéticos. Su construcción duró 
cinco años y costó cinco millones de dólares. 

Unos años más tarde, en 1964, la compañía Digital lanzó al mercado su 
PDP-8, ya equipado con transistores, que se convirtió en el primer miniordenador 
fabricado en serie. Un ordenador de este tipo fue el encargado de la delicada mi- 
sión de guiar la nave Apolo 11, así como su módulo de descenso, en aquel memo- 
rable 21 de julio de 1969 en que, por primera vez, el hombre puso sus pies en 
la Luna. 

Otro de los grandes monstruos que se conservan en el Museo de Boston 
es el ILLIAC IV, utilizado por la Administración Nacional de Aeronáutica y del 
Espacio (NASA) de Estados Unidos. Este equipo estuvo en servicio desde 1975 
hasta 1982, siendo uno de los mayores ordenadores del mundo. Su coste fue de unos 
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100 millones de dólares y necesitaba para funcionar una potencia de medio millón 
de vatios. 

En el campo de la Astronomía, los ordenadores juegan un importante papel 
en la realización de múltiples funciones, tales como: control y manejo de los teles- 
copios y cúpulas; corrección de los efectos de la aberración de la luz; corrección 
de las flexiones mecánicas de los telescopios; control y adquisición de datos obte- 
nidos por diferentes instrumentos, como los fotómetros del visible y del infrarrojo, 
elementos de detección digitales (sistemas Reticon y cámaras CCD), espectrógrafos, 
etcétera; cálculos de órbitas; transformación de los diferentes tipos de tiempos 
utilizados en Astronomía; cálculo, manejo y procesado de los datos obtenidos por 
los diversos instrumentos, etc. 


Cálculo de la órbita y redescubrimiento del 
cometa en su paso de 1986 
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El 21 de octubre de 1982, la Unión Astronómica Internacional daba a conocer 
una noticia importante. El 16 de octubre de 1982, el cometa Halley había sido 
redescubierto por un equipo de astrónomos dirigido por David Jewitt y G. Edward 
Danielson, utilizando una cámara CCD (dispositivo de cargas acopladas) incor- 
porada al gran telescopio de 5 m de Monte Palomar. El Halley, pues, acudía 
puntual a su cita. 

Fue fotografiado como un débil objeto, con una magnitud de 24,3, es decir, 
con un brillo 15 millones de veces menor que el más débil objeto perceptible a 
simple vista. El cometa se encontraba en la constelación de Canis Minor, muy 
cerca de la estrella Procion, a sólo 9 segundos de arco de la posición predicha 
por el astrónomo Donald K. Yeomans, del Jet Propulsion Laboratory (California). 
El cometa estaba situado en ese momento a una distancia del Sol de 11,04 u.a., 
es decir, a 1.652 millones de Km más lejos que Saturno. 

El avance conseguido en este redescubrimiento, respecto al realizado en su an- 
terior presentación de 1910, resultan evidentes si comparamos estos datos con los 
de entonces. 

En su paso anterior, fue redescubierto por Max Wolf, director del Observatorio 
de Heidelberg, el 11 de septiembre de 1909, sol sólo 213 días antes del perihelio, 
que ocurrió el 20 de abril de 1910. En la presentación actual, según hemos dicho 
antes, fue redescubierto cuando aún le faltaban 1.215 días para su perihelio. 

La luminosidad que tenía cuando se fotografió en 1909 era de magnitud 15,5, 
es decir, unas 10.000 veces menos brillante que la estrella más débil visible a simple 
vista. En aquella fecha, el cometa se hallaba desviado 2 minutos de arco de la po- 
sición calculada por Cowell y Crommelin. La distancia a la que se encontraba 
del Sol era de 3,4 u.a., es decir, un poco más allá que el planeta Marte, a 
509 millones de Km. 

El éxito alcanzado en la presentación actual ha sido debido, por una parte, 
a un mejor conocimiento de la naturaleza de los cometas, lo que ha permitido 
introducir en las ecuaciones los efectos producidos por las fuerzas no gravita- 


cionales y, por otro, a los avances tecnológicos logrados por el hombre en los años 
que median entre una y otra presentación. 

El cálculo tan exacto de la órbita del cometa ha sido realizado por el astró- 
nomo Donald K. Yeomans, utilizando un complejo programa de ordenador, en el que 
se tenían en cuenta todas las perturbaciones producidas por los planetas del Sistema 
Solar, así como las fuerzas no gravitacionales del cometa. La fecha calculada 
para su paso por el perihelio es el día 9,3 (U.T.) de febrero de 1986. 

Evidentemente, la complejidad y extensión de las ecuaciones que se han necesitado 
para realizar estos cálculos hubieran hecho imposible su resolución sin la ayuda 
de un potente ordenador. 
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Programa espacial 


para la 
observación del 
cometa Halley 


Los cometas son en la actualidad los cuerpos peor conocidos de nuestro 
Sistema Planetario. Existen todavía muchas dudas sobre su auténtica naturaleza 
y origen, dudas que, muy probablemente, serán disipadas a raíz de la actual visita 
del Halley. 

Un estudio profundo de los cometas puede suponer un gran avance en el cono- 
cimiento sobre el origen y formación de nuestro Sistema Planetario, teniendo en 
cuenta que los cometas provienen de una región muy remota del Sistema Solar, 
donde es muy posible que hayan conservado las propiedades iniciales de la ne- 
bulosa primigenia de donde se formaron el Sol y los planetas. 

Desde el año 1960, los astrónomos se dieron cuenta del gran interés que tenía 
el envío de una sonda espacial hacia un cometa, para realizar un estudio en pro- 
fundidad sobre su naturaleza. 

Había que elegir un cometa que presentara una órbita perfectamente conocida 
y que, además, fuera lo suficientemente activo como para ser capaz de desarrollar 
la cola de plasma. Entre los cometas que cumplían estas condiciones, los astró- 
nomos decidieron seleccionar al cometa Halley en su retorno de 1986. 


Seis sondas espaciales para el estudio del cometa 


Cuatro agencias espaciales están implicadas en los diferentes proyectos que 
actualmente están en marcha para el estudio del Halley. La Intercosmos de la 
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Unión Soviética, la ISAS del Japón, la NASA de Estados Unidos de América y 
la ESA de Europa. 


Proyecto de la NASA 


Desde el año 1961, la NASA comenzó a estudiar las posibilidades de una misión 
espacial para el estudio del Halley y decidió enviar una sonda de órbita retró- 
grada, al igual que la del cometa, con el fin de alcanzarlo en su paso por el perihelio 
y acompañarlo en parte de su recorrido, lo que permitiría estudiarlo con- detalle. 

El coste de este programa era muy elevado y los recortes presupuestarios de 
la NASA impidieron su ejecución. Pero los científicos de la NASA no se resig- 
naron a perder esta oportunidad. Por ello, en el año 1983 decidieron utilizar la 
sonda espacial ISEE-3, lanzada en 1978 para estudiar el cometa Giacobini-Zinner, 
en septiembre de 1985 y el Halley, en marzo de 1986. Para realizar esta misión, 
la sonda ha sido rebautizada con el nombre de ICE (International Cometary 
Explorer). 

Tras unas complicadas maniobras, en enero de 1984 fue modificada su trayec- 
toria con el fin de redirigirla, primero, al encuentro del Giacobini-Zinner y, des- 
pués, hacia el Halley. 

Esta sonda se aproximará al cometa Halley a una distancia de 31 millones 
de Km, el 28 de marzo de 1986. 


Proyecto de Intercosmos 


La agencia Intercosmos de la Academia de Ciencias de la Unión Soviética ha 
enviado al encuentro del Halley dos vehículos espaciales, el Vega 1 y el Vega 2, 
en una misión combinada con el planeta Venus. 

Estos vehículos fueron lanzados el 15 y el 21 de diciembre de 1984 hacia 
Venus, donde tras depositar allí las sondas previstas para el estudio de este planeta 
en junio de 1985, han sido redirigidas hacia el Halley, al que alcanzarán el 6 
y el 9 de marzo, respectivamente. 

El Vega 1 tiene previsto su acercamiento hasta 10.000 Km del núcleo del cometa 
y el Vega 2 a unos 7.000 Km. 


Proyecto de la ISAS 


El Instituto Japonés de Ciencia Astronáutica y del Espacio (ISAS) ha enviado 
también dos sondas al encuentro del Halley. La Planet-A (Sakigake) y la MS-T5 
(Suisei). 

La MS-T5 fue lanzada en enero de 1985 y la Planet-A, en agosto del mismo 
año. La primera pasará a 15 millones de Km del núcleo, atravesando la cola de 
plasma del cometa, y la segunda a 100.000 Km. 


Proyecto de la ESA 


Sin lugar a dudas, el proyecto más ambicioso y atrevido es el de la agencia 
espacial europea, integrada por 11 países, entre los que se encuentra España. 

La ESA ha enviado al encuentro del Halley la sonda Giotto, que fue lanzada 
en julio de 1985. 

Su objetivo es acercarse al núcleo del cometa a sólo 500 Km de distancia, 
con lo que sus cámaras obtendrán imágenes con una resolución de 20 m. 

Muy probablemente, la sonda no sobreviva a la experiencia, debido a los 
fuertes impactos de las partículas de polvo que chocarán contra sus escudos 
protectores a una velocidad de 250.000 Km por hora. Por esta causa, la trans- 
misión de los datos se hará en tiempo real, comenzando a transmitir desde 
4 horas antes de su máximo acercamiento. 

Para llevar a cabo todos estos programas internacionales ha tenido que ser 
realizado un esfuerzo científico sin precedentes en la historia de la Astronomía. 
Pero, no cabe la menor duda, que la empresa merece la pena, pues la humanidad 
va a poder desvelar en 1986 los secretos que, durante tantos siglos, guardan ce- 
losamente estos extraños viajeros del espacio interplanetario. 
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Introducción 
a los programas 


Una vez expuestos los conceptos teóricos fundamentales sobre qué son, dónde 
están y cómo se comportan los cometas, vamos a intentar dar una imagen visual 
que ayude a su mejor comprensión y conocimiento. 

Para ello presentamos cinco programas gráficos y dos numéricos como 
complemento para quienes no se conformen con lo teórico y se decidan por la 
práctica de la observación astronómica. 

Como se vio en la primera parte del libro, los cometas son cuerpos de 
tamaño relativamente pequeño que se mueven alrededor del Sol describiendo 
órbitas de Kepler. El primer programa nos permitirá comprobar cómo son las 
trayectorias de los diversos cuerpos que siguen las leyes de Kepler (pág. 51). En 
general, las órbitas cometarias se caracterizan porque tienen mayor excentricidad 
que las planetarias. 

El segundo programa dejará clara esta diferencia. Sobre las órbitas planetarias 
podremos ver la de cualquier cometa sin más que conocer dos de sus parámetros: 
el semieje mayor y la excentricidad. 

Sin lugar a dudas, el más famoso de todos los cometas es el Halley. El tercer 
programa, “Inclinación de la órbita del Halley”, mostrará su trayectoria desde 
antes de llegar a la órbita de Saturno hasta después del perihelio, y podremos 
comprobar, además, la inclinación de su plano orbital sobre la eclíptica (plano de 
la órbita terrestre). 

El cuarto programa, “Halley”, nos será de gran utilidad para conocer las 
posiciones relativas del Sol, la Tierra y el Halley, dándonos una idea de cómo y 
cuándo podrá ser observado. 
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Aún más importante para la observación será conocer en qué parte de la 
cúpula celeste podremos encontrar cada día al Halley. Para ello, el quinto 
programa, “El Halley sobre el cielo nocturno”, nos dará su posición sobre las 
constelaciones más características que sigue en su trayectoria. 

El primero de los programas numéricos nos dará el valor aproximado de la 
magnitud aparente del Halley, según las previsiones de Morris y Green (la 
magnitud es una medida de la luminosidad de un astro; cada orden de magnitud 
indica que el astro es dos veces y media más luminoso que el siguiente; una noche 
clara nos permitirá ver cuerpos celestes de hasta quinta y hasta sexta magnitud). 

El segundo programa numérico nos dará la altura sobre el horizonte de 
cualquier astro en un determinado instante; su visibilidad será mejor cuanto más 
elevado sobre el horizonte se nos presente. 

En los programas “Halley” y “El Halley sobre el cielo nocturno” se incluyen 
también los listados para el ordenador Spectrum. Los demás programas, dado su 
tamaño pequeño, serán fácilmente traducibles a otros lenguajes BASIC siguien- 
do las instrucciones del apéndice C. 

En cualquier caso, el usuario podrá introducir en los programas todas las 
mejoras que considere oportunas para hacerlos más interesantes o, incluso, 
realizar él mismo sus propios programas; para ello, en el apéndice A encontrará 
algunos de los procedimientos utilizados en este libro aplicados al sistema MSX. 


Programa número 1: Orbitas de Kepler 


Kepler consiguió explicar el movimiento de los planetas mediante tres leyes 
empíricas (v.p.). Un cuerpo que gira alrededor de otro siguiendo estas leyes lo 
hace describiendo una órbita que llamamos “kepleriana”. 

En el programa, el Sol aparece situado en el foco de la órbita elíptica, cuya 
excentricidad será la seleccionada por el usuario. Aparecerá un sprite, que 
representa el astro en estudio, que se moverá dejando tras de sí una estela de 
puntos, que marcarán su trayectoria. 

La velocidad de desplazamiento del astro sigue la segunda ley de Kepler (en 
tiempos iguales se barren áreas iguales). En el perihelio (posición más cercana al 
Sol), la velocidad es mayor que en el afelio (posición más alejada del Sol). Más 
exactamente, la velocidad en cada punto será inversamente proporcional a la 
distancia del Sol y la velocidad angular será inversamente proporcional al 
cuadrado de esa distancia. 

La separación entre los puntos que señalan la trayectoria es proporcional a la 
velocidad del movimiento del sprite. El trazado de la trayectoria nos permitirá 
obtener una imagen gráfica de cómo varía ésta. Para mantener más tiempo la 
imagen podemos utilizar el comando STOP, deteniendo la ejecución del pro- 
grama. 

Si introducimos una excentricidad próxima a cero (por ejemplo, 0.05), obten- 
dremos una órbita característica de un planeta: casi circular y con velocidad casi 
constante en toda su trayectoria. 

Con excentricidades mayores de 0.5 obtendremos una clara diferencia de las 
velocidades del perihelio (rápido) y del afelio (muy lento). Son órbitas de tipo 
cometario. 


Figura 1. Orbita de Kepler 


En el programa se han fijado las dimensiones de la elipse (semieje mayor) y el 
punto en que está centrada. Como el Sol se sitúa en un foco de la elipse, al 
aumentar la excentricidad se situará cada vez en una posición más alejada del 
centro. 
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Detalles técnicos del programa 


El sprite se define en la línea 190 y la trayectoria elíptica se parametriza en las 
líneas 330 y 340, 

El argumento de las funciones seno y coseno será una variable T, que aumenta 
según el contador I, cuyo incremento es la variable J. 

El valor de la J será inversamente proporcional al cuadrado de la distancia al 
foco como corresponde a la velocidad angular: líneas 360 y 370. 

La J será para nosotros el equivalente en un movimiento discreto a la 
velocidad angular en uno continuo, ya que determina el aumento en el ángulo por 
cada ciclo (entre las líneas 400 y 520). 

En la línea 350, para calcular la R (cuadrado de la distancia al foco en cada 
instante), tenemos el parámetro C, que será la distancia entre el centro de la elipse 
y sus focos. 

El otro INPUT es la velocidad, que será un factor constante en el cálculo del 
valor de J y que, por tanto, afectará a la velocidad de ejecución del programa. 

La aproximación del movimiento continuo por uno discreto sólo será buena si 
los incrementos de la variable J son lo suficientemente pequeños. Por eso, las 
entradas de excentricidad y velocidad están limitadas, ya que cualquiera de las 


dos puede hacer que J sea demasiado grande. 


O 170 IF E>.89 THEN 460 
27 x* e 180 SCREEN 2,2 
3” * ORBITA DE KEPLER + 190 SPRITES (1)=CHR$ (H1C) +CHR$ (2H1C) 
4 E +CHRS$ (%H1C) +CHRS (HO) +CHR$ (HO) +CHRS ( 
E 4 HO) +CHR$ (XHO) +FCHR$ (LHO) 

d 200 ” 
A KEYOFF:CLS O 


20 LOCATE 0,3 


220 


? PARAMETROS ORBITALES 


30 PRINT "Este programa muestra una o a 
rbita kepleriana tipica,en funcio 240 ” 

n de su excentricidad." 250 Ai= 90 

40 PRINT 260 B1=?0* ABS (SOR (1-ExE) > 
50 PRINT "La velocidad introducida so 270 J=1 

lo afecta a la duracion de la ejecuci 280 R=3000 

on del programa,no teniendo relaci 290 C=A1*E 

on directacon la velocidad real ni co 300 ” 

n la va- riacion de esta a lo largo 310 CIRCLE (125-C,92),3 
de la or- bita." 320 PAINT (1235—E, 922 

60 PRINT 330 X1=A1*COS(T) 

70 PRINT "La variacion de la velocid 340 Yi=B1x*SIN(T) 


ad del co-meta esta ajustada de forma 


350 


R=(X1+C)+*(X1+C0)+Y1+*Y1 


rigurosa a la situacion real." 360 J=2500*V /R 

80 LOCATE 0,20 370 I=1+J 

90 PRINT "Para comenzar pulse barra e 380 T=1/10 

spaciadora" 390. * 

100 IF STRIG(0)=0 THEN 100 400 PUT SPRITE 1, (X1+105,Y1+91) 
110 CLS 410 


120 INPUT "VELOCIDAD (1-5)"3VW 


420 


PSET (X1+109,Y1+92) 


130 IF V>5 THEN GOSUB 500:GOTO 120 430 IF T>2*3.1416 THEN 520 
140 IF V<1 THEN GOSUB 310:GOTO 120 440 ? 

150 PRINT 450 GOTO 330 

160 INPUT "EXCENTRICIDAD";5E 460 PRINT 
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470 PRINT "EXCENTRICIDAD DEMASIADO AL 530 NEXT 1 


TA" 540 SCREEN O 

480 PRINT 5590 LOCATE 0,10 

490 GOTO 160 560 PRINT "Si quiere volver a comenza 
500 PRINT:PRINT "VELOCIDAD DEMASIADO ry,pulse la”":PRINT TAB(10)3"barra espa 
ALTA" : RETURN ciadora” 

510 PRINT:PRINT "VELOCIDAD DEMASIADO 570 IF STRIG(0)=0 THEN 370 

BAJA" : RETURN 380 CLEAR 

520 FOR I=i TO 1500 590 GOTO 10 
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Programa número 2: Orbita de un cometa 


Los planetas tienen unas órbitas casi circulares. Sólo Mercurio y Plutón tienen 
una excentricidad apreciable (del orden de 0.2). En el caso de los cometas, la 
excentricidad puede llegar a ser casi 1.0, como en el caso del Halley, que es de 
0.967. Por tanto, sus órbitas serán muy excéntricas y de forma muy diferente a la 
de los planetas. 

Este programa nos permitirá comprobar gráficamente esa diferencia. Sobre las 
órbitas de la Tierra, Marte, Júpiter, Saturno, Urano y Neptuno nos permite trazar 
la órbita de cualquier cometa sin más que introducir los valores de su semieje 
mayor y su excentricidad. En la tabla de la figura 2 están los parámetros de los 
cometas más conocidos, y en cualquier anuario astronómico es fácil encontrar los 
de otros muchos. 


y 


Encke 
Grigg-Skjellerup 
Tempel 

D'arrest 
Pons-Winnecke 
Kopff 


Brooks 

Faye 

Comas sola 
Crommelin 

Halley 

C. Herschel-Rigollet 
Tuttle 


Figura 2. Datos de algunos cometas que se pueden utilizar en el programa 


La aproximación por instantes discretos que hemos utilizado en el anterior 
programa sólo es buena cuando la excentricidad no es muy alta. Pero en el 
presente programa, para evitar que la duración de su ejecución sea excesiva y 
permitir que la excentricidad pueda ser alta, hemos preferido que el sprite que 
representa al cometa no siga exactamente la ley de Kepler, sino que se disminuye 
la diferencia de velocidades entre perihelio y afelio. No olvidemos que si lo 
hiciéramos en “tiempo real”, en el caso del Halley nos llevaría setenta y seis años. 


Detalles técnicos del programa 


Los parámetros orbitales se introducen con las sentencias INPUT de las líneas 
180 y 220. 

Entre las líneas 310 y 360 se define el sprite y los valores del resto de 
parámetros del cometa. El Sol y las órbitas planetarias se definen entre las líneas 
420 y 490. 


Figura 3. Orbita del cometa Encke 


Figura 4. Orbita del cometa Tuttle 
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Las posiciones del sprite están dadas por las variables X1 Y1, parametriza- 
das en las líneas 550 y 560 por la variable T. 

La sentencia CLEAR de la línea 800 pone a cero todas las variables para 
poder volver a ejecutar el programa. 


7 II IA IEEE MIMI 
> * 
? * ORBITA DE UN COMETA * 
. + * 
AS 
> 


00 AQNr 


10 CLS 

20 COLOR 15,1,1 

30 PRINT "Este programa muestra la or 
bita de un cometa dentro del sistem 
a solar." 

40 PRINT 

50 PRINT "La velocidad introducida so 
lo afecta a la duracion de la ejecuci 


on del programa,no teniendo relaci 
on directacon la velocidad real ni co 
n la va- riacion de esta a lo largo 

de la or— bita." 

£0 PRINT 


70 PRINT "Tanto las orbitas del comet 

a como lasde la Tierra,Marte, Jupiter 
s¿Saturno, Urano y Neptuno estan repre 

sentadas aescala real.” 

380 PRINT 

90 PRINT "Las velocidades del cometa 

a lo lar- go de su orbita no se corre 

sponden con la autentica diferencia 
de velo- cidades entre el afelio y € 
1 perihe- lio laproximadamente de 1 « 
60)." 

100 LOCATE 0,21 

110 PRINT "Pare comenzar pulse barra 
espaciadora” 

120 IF STRIG(0)=0 THEN 120 

130 CLS 

140 INPUT "VELOCIDAD (0-10)"5V 

150 PRINT 

1650 IF V>10 THEN 710 

170 PRINT 

180 INPUT "SEMIEJE MAYOR (lu.a.) "53 Ac 
190 IF AOX<1 THEN. 180 

200 A1=A0*5 

210 PRINT:PRINT 

220 INPUT "EXCENTRICIDAD ";E 

230 PRINT 

240 1F E>1 THEN PRINT "EXCENTRICIDA! 
DEMASIADO GRANDE":GOTO 210 

250 SCREEN 2,2 

260 ?” 

A 

280 ” SPRITE Y FARAMETROS 

ZOO? IIA MAA AL ME MIE MER AMET 

300 ” 

310 SPRITES (1)=CHR% (%H18) +CHRS (£H18) 
+CHR$ (UH0) +CHR$ (HO) +CHR$ (HO) +CHRS$ (e 
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HO) +CHR$ (%HO) +CHR$ (LHO) 

320 ” 

330 B1=A1+*ABS(SOR (1—E+*E) ) 

340 J=1 

350 C=A1x*E 

360 R=1 

370 ?” 

A AAA 
390 ”? SOL Y ORBITAS PLANETARIAS 
400 ” FEMME DEMI AE MEE DE IE IDE ME DEBIO DIE TIE IE 
410 >” 

420 CIRCLE (57,92),1,11 

430 PAINT  (57,92),11,11 

440 CIRCLE (57,92 ),5,4 ¿REM tierra 
450 CIRCLE (57,92 ),9,4:REM marte 
460 CIRCLE (57,92 ),29,4:REM jupiter 
470 CIRCLE (57,792 ),56,42: REM saturno 
480 CIRCLE (57,92 ),106,4,4.5,1.5: REM 
urano 

490 CIRCLE (57,92 ),167,4:REM neptuno 
300 ” 

A 

520 ? MOVIMIENTO 

A 

540 ” 

550 X1=A1*COS(T) 

560 Y1i=B1*SIN(T) 

570 R=SOR ( (X1+C0)*(X1+C)+(V1x*Y1)) 
580 R=R" (2/3) 

590 J=V/R 

600 I=I+J 

610 T=1/10 

620 ” 

630 1F X1+54+C>255 OR X1+54+C<00R Y1+ 
91>1791 ORY1+91<0OTHEN 550 

640 PUT SPRITE 1, (X1+544+C,Y1+91) 
650 PSET (X1+57+C,Y1+92),4 

660 ” 

670 IF T>2*3.1416 THEN 740 

680 ” 

690 GOTO 550 

700 PRINT "" 

710 PRINT"VELOCIDAD DEMASIADO ALTA" 
720 PRINT"" 

730 GOTO 140 

740 FOR 1l=1 TO 1500 

750 NEXT 1 

760 SCREEN O 

770 LOCATE 1,10 

780 PRINT"Si quiere volver a comenzar 
¿pulse la"¿PRINT TAB(10)3 "barra espac 
iadora” 

790 1F STRIG(0)=0 THEN 790 

800 CLEAR 

810 GOTO 10 


Programa número 3: Inclinación de la órbita 
del Halley 


Salvo para Neptuno y Plutón, el ángulo formado por los planos orbitales de 
los planetas con la eclíptica (plano orbital de la Tierra) es muy pequeño, por lo 
que, en una primera aproximación, podemos decir que sus órbitas son “coplana- 
rias”. Pero no ocurre lo mismo con los cometas. En general, su plano orbital no 
coincidirá con el de los planetas. Así, el plano orbital del cometa Halley forma un 
ángulo de 18 grados con la eclíptica. 

Por medio del presente programa podemos visualizar el efecto de esta 
diferencia de planos desde la perspectiva que se quiera tomar, ya que representa la 
trayectoria del Halley desde que llega a la órbita de Saturno hasta después del 
perihelio. 

Un valor próximo a cero para la inclinación corresponde a una visión “a vista 
de pájaro” del plano de los planetas. La órbita del Halley se proyectará sobre 
ellos, no apreciándose el ángulo de su plano. 

Si introducimos una inclinación cercana a 90 grados, tendremos la imagen 
más demostrativa del ángulo entre ambos planos orbitales. A 90 grados, al 
aparecer la eclíptica de perfil, se visualizará como una línea hacia la cual se 
dirige el cometa en una trayectoria en evidente ángulo de unos 18 grados. 

En todas las inclinaciones, en el momento de que el plano de la órbita 
cometaria atraviesa la eclíptica, hacemos un cambio de color de la trayectoria, con 
lo que puede fácilmente trazarse la línea de corte entre ambos planos. 

Para poder mantener la posición correcta entre velocidades y posiciones del 
Halley con los demás planetas, ha sido necesario que la velocidad del cometa en 
posiciones alejadas del perihelio fuera muy pequeña. 


Figura 5. Inclinación de la órbita del Halley. Inclinación 60" 
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Figura 6. Inclinación de la órbita del Halley. Inclinación 90% 


Detalles técnicos del programa 


Las órbitas se dibujan por medio de las sentencias CIRCLE de las líneas 490 


y 520. 


La excentricidad de las órbitas será A, que se obtiene en las líneas 160 y 170 


a partir del valor introducido para la inclinación en la línea 130. Este valor de 
A también afectará a la parametrización de la trayectoria del Halley en la lí- 
nea 770. Tenemos que aplicar una transformación del sistema de coordenadas que 
definen la posición del cometa para corregir la inclinación real de su órbita, líneas 
790 y 800, donde pasamos de las variables X0, YO a X1, Y1, que serán las 
utilizadas en la línea 910 para calcular el punto Z1, Z2, en el que debe colo- 


carse el sprite, línea 920. 


Para que coincida el centro del sprite con la trayectoria, lo pondremos en 
(Z1-3, Z2-2). Recordemos que la sentencia PUT SPRITE (X, Y) nos pone en X, Y 
el extremo superior izquierdo del sprite. 


d, A RS 
2 * 
3 ?* ORBITA DEL COMETA HALLEY  * 
4 >% * 
O A 
6 > 

10 CLS 


20 COLOR 15,1,1 

30 LOCATE 0,1 

40 PRINT "Este programa muestra la or 
bita del Halley dentro del sistema s 
olar."” 

50 PRINT 

60 PRINT "Las velocidades,posiciones, 
y las or— bitas de Saturno, Jupiter , 
Marte y Tierra,asi como los demas p 
rametros estan realizados a escala r 
eal,siendo” 

70 PRINT "por tanto, todas las posici 
ones rela—-tivas de los planetas y el 

cometa lascorrectas en cada momento.” 
80 LOCATE 0,20 
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90 PRINT"Para comenzar pulse barra es 
paciadora” 

100 IF STRIG(0)=0 THEN 100 

110 CLS 

120 V=10 

130 INPUT "Grados de inclinacion”;A0 

140 IF A0>90 THEN CLS:GOTO 130 

150 E=.96728 

160 AO=40%*3.14159/180 

170 A=COS(A0) i 

180 ” 

A 

200 ” DEFINICION DE SPRITES 

A 

220 ” 

230 SCREEN 2,2 

240 SPRITES (1)=CHR$ (xH18)+CHR$ (2H18) 
+CHR$ (YHO) +CHR$ (%HO) +CHR$ (HO) +CHRS (Ze 
HO) +CHR$ (HO) +CHR$ (HO) 

250 SPRITES (2) =CHR$ (2%H1C) +CHR$ (XH1C) 
+CHR$ (%H1C) +CHR$ (HO) +CHR$ (HO) +FCHRS ( 
HO) +CHR$ (%HO) +CHR$ (HO) 


260 SPRITES (3)=CHR$ (HO) +CHR$ (H3C) + 
CARS (HSA) +CHR$ (RH3C) +CHR$ (LHO) +CHR$ ( 
HO) +CHR$ (HO) +CHR$ (HO) 
270 ” 
A 
z20 > PARAMETROS DEL HALLEY 
O 
310 ” 
320 Al=177 
330 Bl1=177*ABS(SOR (1-E*E) ) 
340 C=A1x*E 
350 HX=COS (18*(1-A)%3.14/180) :¿ HY=SIN( 
-18* (1-4)+*3.14/180) 
360 R=10 
3/0 ” 
A 
390 ?” EL SOL Y.LA ORBITAS 
7 MARIA III AAA 


420 ”El Sol 

430 ? 

440 CIRCLE (120,95),2,11 

450 PAINT  (120,95),11,11 
> 


470 ”*Las orbitas 
480 ?” 
490 CIRCLE (120,95),11,4,,»A:REM TIERR 


500 CIRCLE (120,95),16,4,,,AM: REM MARTE 
510 CIRCLE (120,95),51,4,,,AM:REM JUPIT 


520 CIRCLE (120,95),794,4,,,A2REM SATUR 


330 ” 

A 
300 ? CONTADORES 
A 
570 ” 

580 J=V/R 

590 I=I+J 


600 T=1.05+1 
610 L=L+1.6E-03 


FEDEAHEMEHEIO IE MED METE DE MIEDO EE DEE IE IE METE MEE IEEE dE 


ECUACIONES DEL MOVIMIENTO 
VEMM IE MMM DE MMM IE MEM IEA IRE 


o 
+ 
[e] 
. “uo os 


670 X2=10*C0S (-2. 47+76*L) 2: Y2=-10*SIN ( 
—-2.47+76*L) 

680 X3=13*C0S (—. 6+490. 44*L) : Y3=-15*SIN 
(—.6+40. 44+*L) 

690 X4=51*C0S (2.37+6.4*L) : Y4=-5S1%4SIN ( 
2.37+6. 4*L) 

700 X3=94%*C0S (2. 58+L) : YS=-94*SIN (2.58 
*L) 


710 ” 
A AA 
730 ” CAMBIO DEL SISTEMA 


DE COORDENADAS Y 
TRAYECTORIA DEL HALLEY 
A 
Za * 
760 XO=A1+*C0S(T) 
770 YO=B1i*(A+.05S)*SIN(T) 
780 Yó6=B1*SIN(T) 
790 X1=XO*HX+YO*HY 
800 Y1=-XO*HY+YO*HX 
310 R=(X0+C)+*(X0+C)+(Y56*Y6) 
820 ” 
7 HEM RIR RIIIE 
840 >” UBICACION DE LOS SPRITES 
7 MMMM MIMI III MEM IEAAM 
860 ” 
870 PUT SPRITE 2, (X2+117,Y2*A+94),,1 
880 PUT SPRITE 3, (X3+117,Y3*AM+94),,1 
890 PUT SPRITE 4, (X4+117,Y4x*A+94),,2 
900 PUT SPRITE 5, (X5+118,Y5*A+93),,3 
910 Z1=X1+288-Cx* (1-HX) 2 Z2=Y1+94-Cx*HY 
920 PUT SPRITE 1, (Z1-3,Z2-2) 
930 IF Z1>120 THEN K=6 ELSE K=3 
940 PSET (Z1,22),K 
950 IF T>3.93THEN 970 
960 GOTO 580 
970 FOR I=1 TO 1500 
980 NEXT I 
990 SCREEN O 
1000 LOCATE 1,10 
1010 PRINT"Si quiere volver a comenza 
ry,pulse la":PRINT TAB(10)5 "barra espa 
ciadora" 
1020 IF STRIG(0)=0 THEN 1020 
1030 CLEAR 
1040 GOTO 10 
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Programa número 4: “Halley” 


Las condiciones de visibilidad del cometa Halley desde la Tierra dependerán 
tanto de las posiciones relativas Tierra-Sol-cometa, como del tamaño y orienta- 
ción de la cola de éste. 

En este programa podremos ver la evolución de estos factores desde julio de 
1985, cuando desde la Tierra se lanza la sonda Giotto y el cometa aún no ha 
alcanzado la órbita de Marte, hasta que después de haber pasado el perihelio 
vuelva a alejarse de nosotros. 

Cuando el Halley aún no haya llegado a la órbita de Marte, su cola será 
pequeña y, dada la distancia que lo separa de nosotros, apenas será visible. 

Observando con detenimiento la pantalla podemos, en todo momento, hacer 
un razonamiento similar que nos permitirá hacer una estimación del aspecto que 
presentará el Halley en cada situación. Por otro lado, pulsando la barra espacia- 
dora, obtendremos en pantalla una tabla con los datos más interesantes en cada 
parte de su trayectoria. 

Por ejemplo, es fácil ver que la distancia Tierra-cometa pasa por un mínimo 
antes del perihelio y, por otro, después de él. Tendremos entonces dos... mínimos 
relativos en su magnitud. También resulta evidente que el cometa no será visible 
en el perihelio, pues estará en línea con el Sol. 

Fijándonos un poco más, igualmente será posible conocer las horas aproxima- 
das en que será observable el cometa, considerando el ángulo Tierra-Sol-Halley. 
Antes del perihelio, el cometa estará al Este del Sol y, por tanto, será visible un 
tiempo después del anochecer (mayor cuanto mayor sea el ángulo Tierra-Sol- 
Halley). Pasado el perihelio, el Halley estará al Oeste del Sol y, por tanto, saldrá 
un tiempo antes de éste y sólo será observable un poco antes del amanecer. Este 
tiempo aumentará al aumentar el ángulo entre los tres astros. 

Los datos que aparecerán en pantalla al pulsar la barra espaciadora son los 
previstos actualmente. Si hay cambios, pueden ser corregidos en el listado (líneas 
2370-2500). 


Detalles técnicos del programa 


La sentencia CLEAR de la línea 20 nos aumenta el área de memoria dedicada 
a caracteres hasta 500 bytes; de esta forma, se asegura el poder ampliar la tabla 
de datos con los que resulten de la observación directa. 

Fijados los valores de algunas constantes (líneas 40-70), y tras los textos de 
presentación, se definen los sprites de los planetas y del cometa entre las líneas 
330 y 890. Para ilustrar de la forma más real posible los distintos aspectos que 
presentará el Halley a lo largo de su trayectoria, se definen trece sprites con 
diferentes orientaciones y tamaños de cola. La figura 7 nos muestra algunos 
diseños de estos sprites. 

Las órbitas de los planetas y del Halley se trazan mediante las sentencias 
CIRCLE de las líneas 1130 y 1170, con un factor común de excentricidad, E 
(definido en línea 40), que mos modula la inclinación del ángulo de visión. 


y be de 
ral alan 


Figura 7. Definición de los sprites que representan al cometa Halley 


Para mover los sprites haremos un bucle entre las líneas 1230 y 2100, que nos 
aumentará el valor de los contadores J (línea 1230) y T (línea 1250). 

En la línea 1340 se abre un fichero para poder escribir en la pantalla de 
gráficos. 

A partir de la línea 1410, y hasta la 1600, se definen las ecuaciones del 
movimiento de los planetas, cometa y de la sonda Giotto. 

Las velocidades de cada planeta y de la sonda están moduladas por los 
denominadores y la posición inicial por las constantes, que aparecen en la 
expresión de las variables L. 

Se ha dividido la trayectoria del Halley en trece tramos, en cada uno de los 
cuales es representado por un sprite y se le asigna un valor de Q. Esta variable 
modula la velocidad del cometa en los diferentes puntos de su trayectoria. 

Esta aproximación discontinua será válida, ya que cada tramo es una parte 
muy pequeña de la órbita. 

Los valores de la variable R seleccionan el sprite del Halley para cada posición 
mediante las sentencias IF-THEN de las líneas 1890 a 2080. 

Las variables X5, Y5 definen el movimiento de la sonda Giotto. El radio de su 
trayectoria disminuye (línea 1260) desde un valor inicial, igual al de la Tierra, 
hasta encontrarse con el Halley, cerca de la órbita de Venus. Posiblemente la 
sonda se destruirá al contactar con el cometa (cuando se cumplan las condiciones 
de la línea 1770) y, por tanto, desaparecerá de la pantalla. 

La línea 2090 nos permite, pulsando la barra espaciadora, acceder en cada 
momento a la tabla de datos definida entre las líneas 2370 y 2480. 

Las expresiones de las variables de caracteres X1$-X5$ dependerán del valor 
de R (que indica el punto de la trayectoria), según se define entre las líneas 2180 
y 2300. 
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Figura 8. El Halley antes de llegar a la órbita de Marte 


Figura 9. El Halley en el perihelio 


La función STRIG(0) nos permite volver a la pantalla de gráficos cuando se 
pulsa la barra espaciadora (STRIG(0)= — 1). 

Modificando los valores de las constantes que intervienen en las ecuaciones del 
movimiento (contadores, excentricidad, etc.) se obtendrán diferentes visiones. Pero 
puede que entonces no haya una correspondencia exacta con la realidad. 


1 REM III 
2 REM + * 
3 REM * HALLEY + 
4 REM + Je 
S REM 41330732 
6 REM 

10 CLS:KEYOFF 

20 CLEAR 500 

30 COLOR 15,1,1 


40 E=.6 

SO J=-3 

60 V(0)=2 

70 Q=5*V (0) 

80 ” 

DO? MIAMI MM AMIA 
100 ? TEXTO DE PRESENTACION 


A 

120 ” 

130 SCREEN O 

140 PRINT "Este programa ha sido pens 
ado para representar lo mas fielmen 
te posible el encuentro del cometa Ha 
lley con laTierra en 1986." 

150 PRINT 

150 PRINT "Todas las posiciones relat 


ivas de plan=tas, cometa y Sol,asi 
como sus velocidades, son reales.” 
170 PRINT 


180 PRINT "En cualquier momento, pulsa 

ndo la barra espaciadora, obtendr 
a los para-metros y magnitudes import 
antes del cometa en ese instante y e 

n esa posi-cion." 

190 LOCATE 0,20 

200 PRINT TAB(7)3"PULSE BARRA ESPACIA 
DORA" 

210 IF STRIG(0)=0 THEN 210 

220 CLS 

230 PRINT "El sistema solar esta repr 
esentado enuna perspectiva de unos 45 
grados.” 

240 PRINT 

2530 PRINT "El programa abarca un peri 
odo de docemeses, comenzando desde el 


15 de agosto de 17985 hasta el 15 
de agosto de 1984." 
260 PRINT 


270 PRINT "Recuerde que en cualquier 
momento puede obtener los datos re 
ferentes almovimiento en el instante 
en que pul-se la barra espaciadora, p 
udiendo continuar “cuando lo desee 
desde ese mismo punto." 


280 LOCATE 0,20 

290 PRINT "PARA COMENZAR PULSE BARRA 
ESPACIADORA" 

300 IF STRIG(0)=0 THEN 300 


A 
330 ”? DEFINICION DE SPRITES 
A 
350 ” 

360 GOSUB 2550 

370 CLS:SCREEN 2,2 

380 REM posiciones del Halley 


400 A1$=CHR$ (HO) +CHR$ (HO) +CHR$ (LHO) 
+CHR$ (KUHO) +CHR$ (HO) +CHR$ (HO) +CHRS (2 
HO) +CHR$ (HO) 

410 C1$=CHR$ (2H0) +CHR$ (HO) +CHR$ (HO) 
+CHR$ (KHO) +CHR$ (HO) +CHR$ (KHO) +CHR$ (Le 
H38) +CHR$ (%HCO) 

420 SPRITES (13)=A1$+41$+C01$+A15$ 

430 C2$=CHR$ (HO) +CHR$ (HO) +CHR$ (HO) 
+CHR$ (%H8) +CHR$ (2H1 4) +CHR$ (%H28) +CHR$ 
(%£H60 ) +CHR$ (%H80) 

440 SPRITES (12)=41$+A1$+C2$+1$ 

450 C3$=CHR$ (HO) +CHR$ (2:H20) +CHR$ (2H1 
0) +CHR$ (%H28) +CHR$ (LH5O) +CHR$ (%H60) +C 
HR$ (XHCO) +CHR$ (%H8O) 

460 SPRITES$(11)=A1$+A1$+C3$+15$ 

470 A4$=CHR$% (UHO) +CHR$ (YH 1) +CHRS (2%:H2) 
+CHR$ (%H1 ) +CHR$ (HO) +CHR$ (2H 1 ) +CHR$ (8 
HO) +CHR$ (2%H0) 

480 C4$=CHR$ (HO) +CHR$ (4H40) +CHR$ (HA 
0) +CHR$ (H40) +CHR$ (H80) +CHR$ (3:H40) +C 
HR$ (%H80) +CHR$ (%H80) 

490 SPRITES$ (10)=A4$+41$+C4$+AM1$ 

500 AS$=CHR$ (%:H28) +CHR$ (2H 1 4) +CHR$ (LH 
28) +CHR$ (%H1 4) +CHR$ (HA) +CHR$ (%H6) +CH 
R$ (%H3) +CHR$ (2H 1) 

510 SPRITES (9) =A5$+41$+41$+A1$ 

520 A6$=CHR$ (HO) +CHR$ (HO) +CHR$ (LHAO 
) +¿CHR$ (%H50) +CHR$ (HAS) +CHR$ (2:H54) +CH 
R$ (X%H6) +CHR$ (LH3) 

530 SPRITES (8) =A6$+41$+A51$+A1 $ 

540 A7$=CHR$ (HO) +CHR$ (HO) +CHR$ (HO) 
+CHR$ (HO) +CHR$ (HO) +CHR$ (HO) +CHR$ (2 
H80) +CHR$ (2H54) 

550 B7$=CHR$ (HAB) +CHR$ (4H54) +CHR$ (2H 
390) +CHR$ (HO) +CHR$ (3HO) +CHR$ (HO) +CHR 
$ (XHO) +CHR$ (%HO) 

560 SPRITES (7) =47$+B7$+41$+11$ 

570 B8$=CHR$ (%H3) +CHR$ (%H6) +CHR$ (32H54 
) +CHR$ (KHA8 ) +CHRS$ (%H50) +CHR$ (LXHAO) +CH 
R$ (HO) +CHR$ (HO ) 
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580 SPRITES (6) =A1$+B8$+41$+815$ 

590 BO9$=CHR$ (XxH1 ) +CHR$ (HI) +CHR$S (UHG) 
+CHR$ (HA) +CHR$ (3H1 4) +¿CHR$ (4H28) +CHR$ 
(2414) +CHR$ (2H28) 

600 SPRITES (5) =A1$+B9$+415$+A15$ 

610 B10$=CHR$ (HO) +CHR$ (HO) +CHR$ (2%H1 
) +CHR$ (LHO) +CHRS (%H1)+CHR$ (%H2) +CHR$ ( 
2% H1)+CHR$ (HO) 

620 D10$=CHR$ (%H80) +CHR$ (%H80 ) +EHR$ (Ze 
H40) +CHR$ (£xH80) +CHR$ (£H40) +CHRS (8HAO) 
+CHR$ (%H40) +CHR$ (2%HO) 

630 SPRITES (4) =A1$+B10$+41$+D10$ 

640 D11$=CHR$ (%H80 ) +CHR$ (UHCO) +CHR$ (2: 
H60) +CHR$ (%H50) +CHR$ (4H28 ) +CHR$ (2%:H1 4) 
+CHR$ (%H20) +CHR$ (HO) 

630 SPRITES (3)=41$+41$+415$+D11$ 

660 D12$=CHR$ (%H80) +CHR$ (%xH60) +CHR$ (% 
H28) +CHR$ (4H1 4) +CHR$ (%H9) +CHR$ (HO) +C 
HR$ (HO) +CHR$ (RLHO) 

670 SPRITES (2) =A1$+41$+41$+D125 

680 D13$=CHR$ (RKHCO) +CHR$ (2H38 ) +CHR$ (Re 
HO) +CHR$ (HO) +CHR$ (HO) +CHR$ (2H0) +CHR 
$ (%HO) +CHR+$ (HO) 

690 SPRITES (1)=A1$+A415+41$+D13$ 

700 ? 

E E 

720 *SPRITES DE LOS PLANETAS 

A 

740 ? 

750 REM planeta Mercurio 

760 ? 

770 SPRITES (31) =CHR$ (HO) +CHR$ (XHO)+C 
HR$ (%H10) +CHR$ (U%H38) +CHR$ (%H10) +CHRS ( 
HO) +CHR$ (XHO) +CHR$ (HO) 

780 ” 

790 REM planetas Venus y Tierra 

800 ” 

810 SPRITES (32) =CHR$ (HO) +CHR$ (£H38) + 
CHR$ (%H7C) +CHR$ (UH7C) +CHR$ (HTC) +CHRS 
(4H38) +CHR$ (HO) +CHR$ (HO) 

820 ” 

830 REM planeta Marte 

840 ? 

850 SPRITES (33) =CHR$ (xHO) +CHR$ (HO) +C 
AR$ (2H19) +CHR$ (LH3C) +CHR$ (HC) +CHRE £ 
2H18)+CHR$ (HO) +CHR$ (HO) 

860 ” 

560 SPRITES (7) =A7$+B7$+41$+A15$ 

570 B8$=CHR$ (RH3) +CHR$ (2H) +CHRS (2H54 
) +CHR$ (HAB) +CHR$ (LHSO) + CHAS (LHAO) +CH 
R$ (HO) +CHR$ (HO) 

580 SPRITES (6) =41$+B8$+41 $+41$ 

590 B9$=CHR$ (UH 1) +CHRS (LH3) +CHR$ (UHO) 
+CHR5$ (HA) +CHR$ (2H 14) +CHR$ (:H28) +CHR$ 
(£H14) +CHR$ (2H28) 

600 SPRITES (5) =41$+B9$+41%$+41$ 

610 B10$=CHR$ (%H0) +CHR$ (HO) +CHR$ (2H 1 
) +CHRS (HO) +CHRS (LH 1) +CHRS (%H2) +CHR$ ( 
24H 1) +CHR$ (HO) 

620 D10$=CHR$ (2H80) +CHR$ (%H30) +CHR$ (Le 
H40) +CHR$ (%H80) +CHR$ (2%H40) +CHRS$ (UHADO) 
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+CHR$ (3H40) +CHRS$ (XLHO) 

630 SPRITES (4) =A41$+B10$+41$+D10% 

640 D11$=CHR$ (%H80) +CHR$ (UHCO) +CHRS (% 
H60) +CHR$ (ZHSO) +CHR$ (£:H28) +CHR$ (XH1 4) 
+CHR$ (%H20) +CHR$ (HO) 

650 SPRITES (3)=41$+A1$+41$+D115$ 

660 D12$=CHR$ (£H80) +CHR$ (%H40) +CHR$ (2 
H28) +CHR$ (2H1 4) +CHR$ (%H8) +CHRS (HO) + 
HR$ (2H0) +CHR$ (HO) 

670 SPRITES (2) =A1$+A1$+01$+D125$ 

680 D13$=CHR$ (%HCO) +CHR+$ (%<H38) +CHRS$ (L 
HO) +CHR$ (HO) +CHR$ (HO) +CHRS (HO) +EHR 
$ (KHO) +CHR$ (HO ) 

690 SPRITES (1) =A1$+A41$+41$+D13$ 


A 
720 ”*SPRITES DE LOS PLANETAS 
A 
740 ? 

7530 REM planeta Mercurio 


770 SPRITES (31) =CHR$ (HO) +CHR$ (HO) +€ 
HR$ (£H10) +CHR$ (£H38) +CHR$ (%H10) +CHR$ ( 
£HO) +CHR$ (HO) +CHR$ (HO) 

780 ” 

770 REM planetas Venus y Tierra 

800 ” 

810 SPRITES (32) =CHR*$ (HO) +CHR$ (2:H38) + 
CHR$ (£H7C) +CHR$ (%H7C) +CHR$ (UHTC) +CHR$ 
(2H38) +CHRS$ (HO) +CHRS (HO) 

820 ” 

830 REM planeta Marte 

840 ?” 

850 SPRITES (33) =CHR$ (HO) +CHR$ (LHO.) +C 
HRS (2H18) +CHR$ (£H3C) +CHR$ (%H3C) +CHR$S ( 
2H18)+CHR$ (2H0) +CHR$ (BHO) 

860 >” 

870 REM Giotto 

g80o ” 

890 SPRITES (15) =CHRS$ (2H18) +CHR$ (XxH18) 
+CHR$ (LHO) +FCHR$ (LHO) +CHR$ (RHO) +CHR$ (2 
HO) +CHR$ (HO) +CHR$ (2:H0 ) 

900 ?” 

O 

920 ”? RADIOS DE LAS ORBITAS 

A 

940 ?” 

950 At1)=25 ¿REM Mercurio 

960 A(2)=47 ¿REM Venus 

970 A(3)=46 ¿REM Tierra 

980 A(4)=100:REM Marte 

990 A(0)=1180 :B(0)=299: REM Halley 
1000 A(5)=62 

1010 ?” 

A 

1030 ”? DIBUJO DE LAS ORBITAS 

O 

1050 ” 

1060 REM el Sol 


1080 CIRCLE (115,75),4,11 


1090 PAINT (115,95),11 

1100 >” 

1110 REM orbitas planetarias 

1120 >” 

1130 CIRCLE (119,95),A(1),13,,,E 

1140 CIRCLE (114,95),A4(2),13,,, E 

1150 CIRCLE (115,95),4A(3),13,,,E 

1160 CIRCLE (123,95),A(4),13,,,€ 

1170 CIRCLE (1257,95),A(0),13,,,E*B(0) 

/A(O) 

1180 >” 

AAA 

1200 ”? CONTADOR Y VELOCIDAD 

1210 ? IMM II III IR IE 

1220 ” 

1230 J=J+5 

1240 I=J+*.122x*V (0) 

1250 T=T+Q 

1260 A(S)=A(5)+*(1+710)/(1+711) 

1270 ” 

A 

1290 ? MOVIMIENTO PLANETARIO 

O 

1310 ” 

1320 PRESET (3,160) 

1330 IF J>1 THEN 1340 

1340 OPEN "GRP: "FOR GUTPUT AS +1 

1350 GOSUB 2550:PRINT +1," La sonda 

Giotto es lanzada des 

de la Tierra el 3 

de Julio de 1985.”: GOSUB 

2550 

1360 IFZ=0 OR K=1 THEN 1410 

1370 LINE (4,160)-(160,190),1,BF 
1380 PSET (5,160) 

1390 GOSUB 2550:PRINT $1," La sonda G 


iotto se encuentra co 

n el Halley el 13 de 
Marzo de 1986. ": GOSUB 

2550 

1400 K=1 


1410 REM Marte 

1420 L(4)=-.1351+1/68.698 

1430 X(4)=A (4) *COS(L (4) ) : Y (4) =-A (4) 5 
IN(L (4)) 

1440 REM Venus 

1450 L(2)=-1.6725+1/22.497 

1460 X(2)=A (2) *C0S (L (2)) : Y (2) =-A(2)x*S 
IN (L (2)) 

1470 REM Tierra y Giotto 

1480 G=3.3278+1/28.5 

1490 L(3)=3.3278+1/36.526 

1500 X(5)=4(5)*C0S (6) : Y (5) =-A (5) *SIN( 
G) 

1510 X(3)=4(3)x*C05S (L (3)) 2: Y (3) =-A (3)+S 
IN (L (3)) 

1520 REM Mercurio 

1530 L(1)=2.6417+1/8.797 

1540 X(1)=A(1)x*CO0S(L(1)):2Y1(1)=-A(1)%S5 
IN (L (1)) 

1550 ” 


1560 >” 
1570 REM movimiento del HALLEY 
1580 ?” 
1590 L(0)=2.5+T*V(0)/2000 
16500 X(0)=A (0) *C0S (L(0)): Y(0)=B(0)%SI 
N(L (0)) 
1610 ” 
A 
1530 ”? PONER LOS SPRITES 
a E 
16550 ” 
15660 REM mercurio 
16570 ” 
1580 PUT SPRITE 1, (X(1)+116,Ex*Y(1)+91 
),14,31 
1590 >” 
1700 REM Venus 
1710 >” 
1720 PUT SPRITE 2, (X(2)+111,ExY(2)+91 
de e : 
1730 1F Z=1 THEN 1790 
1740 ” 
1750 REM Tierra y Giotto 
1760 ” 
1770 IF X(5)+101<H AND Y(0)<0 THEN X( 
5)=0: Y(5)=0:Z=1 
1780 PUT SPRITE 15, (X(5)+109,Ex*Y (5)+9 
4),15,15 
1790 PUT SPRITE 3, (X(3)+112,Ex*Y (3)+91 
y E 
1800 ?” 
1810 REM Marte 
1820 >” 
1830 PUT SPRITE 4, (X (4)+120,ExY (4)+91 
1.9133 
1940 ” 
A 
1860 >” SELECCIONAR EL SPRITE 
QUE REPRESENTA AL HALLEY 
A 
1880 ” 
1890 IF R=0 THEN R=1 
1900 H=X (0)+1:249 
1910 IF Y(0)>0 THEN 2000 
1920 IF H>200 AND H<256 THEN R=13:0= 


1930 IF H>151 AND H<200 THEN R=12:0=4 
1940 IF H>114 AND H<151 THEN R=11:0Q 


1950 IF H>97 AND H<114 THEN R=10:0=6 


1960 IF H>82 AND H<97 THEN R=9:Q=9.5 
1970 IF H>71 AND H<82 THEN R=8:0Q=13. 


1980 IF H>255 THEN 2610 

1990 GOTO 2080 

2000 IF  H<71 THEN R=7:0=17.5 

2010 IF H>71 AND H<87 THEN R=6:0=13.5 
2020 IF H>87 AND H<97 THEN R=5:0=9.5 

2030 IF H>97 AND H< 114 THEN R=4:0=6. 
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2040 IF H>114AND H<155 THEN R=3:0=5 

2050 IF H>155AND H<207 THEN R=2 :Q=4 

2060 1F H>207AND H<256 THEN R=1: 0=3. 

5 

2070 IF H>255 THEN 1230 

2080 PUT SPRITE 0, (X10)+1249,Ex*Y(0)+8 

7),15,R 

2090 IF STRIG(0)=-1 THEN 2160 

2100 GOTO 1230 

2110 ” 

BAZO ? HAMMER 

2130 > TABLA DE DATOS 

A 

2150 ” 

2160 CLS 

2170 SCREEN O 

2180 IF R=1 THEN X1$="15-8-85 y 30-9- 

85":X2$="13.4 y 11.2":X3$="6h y óh 12 

m":X4$="19q 7m y 20q":X5$="400" 

2190 IF R=2 THEN X1$="30-9-895 y 24-10 

-85":X2$="11.2 y 9.5":X3$="6h 12m y 5 

h 46m":X4$="20qg y 21g 12m":X5$="280" 

2200 1F R=3 THEN X1$="24-10-85 y 15-1 

2-85"2X2$="9.5 y 6.2":X3$="5h 46m Y 2 

3h 16m":X4$="21g 12m Y 3qg 37m":X5$="1 
10" 

2210 IF R=4 THEN X1$="15-12-85 y 7-1- 
86":X2$="6.2 y 5.6":X3$="23h im y 22h 
2m":X4$="3q 37m y -3q 53m":X5$="150" 
2220 1F R=S THEN X1$="7-1-86 21-1-86" 
:X26="5.6 y 4.8":X38="22h 2m y 21h 37 
m":X4$="-3qg 53m y -69 28m": X5$="200" 
2230 IF R=6 THENX1$="21-1-86 y 4-2-86 
"2X2$="4.8 y 4":X3$="21h 37m y 21h 12 
m":X4$="-6g 28m y -9q 15m":X5$="225" 
2240 IF R=7 THEN X1$="4-2-86 y 15-2-8 
6":X2$="4 y 4.1":X3$="21h 12m y 20h 5 
im":X4$="-9g 15m y -159 59m": X5$="223 


2250 IF R=8 THEN X1$="15-2-86 y 1-3-8 
6":X24$="4.1 y 4.4":X3$="20h Sim y 20h 
26m":X4$="-11g 59m y —16g 19m":X5$=" 

200" 

2260 IF R=9 THEN X1$="1-3-B86 y 21-3-8 
6":X2$="4,4 y 4.4":X3$="20h 26m y 19h 
38m":Xx4$="-1b6gq 19m y -26g 43m": X5$4$=" 
150" 

2270 1F R=10 THEN X1$="21-3-86 y 9-4- 
86":X2$="4.4 y 3.9":X3$="19h 38m y 15 
h 48m":X4$="-26g9 42m y -47g 23m":X5$= 
0) 80 '” 

2280 IF R=11 THEN X1$="9-4-86 y 8-5-8 

6":X2$="3.9 y 6.8":X3$="15h 48m y 10h 
38m":X4$="-47g im y —13g 26m": X3$="1 

05 ” 

2290 IF R=12 THEN X1$="8-5-86 y 15-6- 

86":Xx2$="65.8 Y 9.4": X3$="10h 38m 10h 
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26m”":X4$="-13q 26m Y -5g 21m":X5$="2 
25" 
2300 IF R=13 THEN X1$="15-6-86 y 15-8 
-86":X2$="9,4 y 11.6":X3$="10h 27m y 
10h 59m":X4$="-5g 21m y —6g 39m": X5$= 
n 400 . 
2310 ” 
ZIZO ? MIA AMI MI IA MITIN IE IE RIEIAMA 
2330 ” PRESENTACION DE DATOS 
A A 
2350 ” 
2360 GOSUB 2550 
2370 PRINT "Estamos entre "3X1% 
2380 LOCATE 0,3 

2390 PRINT "MAGNITUD APROXIMADA: Entr 
e "3X26 
2400 LOCATE 0,6 

2410 PRINT "COORDENADAS CELESTES: " 
2420 PRINT 
2430 PRINT ” A.R.:zEntre "3X3$8 

2440 PRINT " D.:Entre "3X4%$ 
2450 LOCATEO, 12 
2460 PRINT "DISTANCIA APROXIMADA A LA 
TIERRA (en millones de km.) "3X5%6 
2470 LOCATE 0,20 

2480 PRINT "PARA CONTINUAR PULSE BARR 
A SEPARADORA" 
2490 IF STRIG(0)=0 THEN GOTO 2490 
2500 GOTO 360 


A 


2530 ? RUTINA DE SONIDO 
O 


2350 PLAY "T10051M5000L46407CCCC":RETU 


2560 ” 
a AAA 


2580 ” TEXTO FINAL 

A 
2600 ” 

2610 CLS 


2620 SCREEN O 

2630 PRINT "Este ha sido el recorrido 
del cometa Halley en_su encuentro co 
n la Tierra correspondiente a 1986." 
2640 PRINT 

2650 PRINT "Todos los parametros y di 
mensiones delos graficos estan repres 
entados es- trictamente en escala,exc 
epto los ra-dios planetarios y del So 
ljasi como la longitud de la cola de 

1 cometa." 

2660 LOCATE 0,19 

2670 PRINT "Para volver a empezar,”":P 
RINT "pulsar la barra espaciadora." 
2680 GOSUB 2550: FOR T=1 TO 1000:NEXT 
T:1F STRIG(0)=0 THEN 2680 

2690 GOTO 10 


El programa “Halley” (versión Spectrum) 


El programa Halley versión Spectrum engloba las dos partes de la versión 
MSX en una: a medida que transcurre la órbita del cometa, los datos de sus 
características van saliendo en pantalla. Esta versión presenta, además, una 
notoria diferencia, ya que, al carecer el Spectrum de sprites, las figuras de planetas, 
Halley y la sonda Giotto, que desde la Tierra sale al encuentro de éste, se 
consiguen con puntos (sentencias Plot). 

Al principio, tras el gráfico y textos de presentación se dibujan las órbitas 
planetarias y del cometa, que son las elipses definidas entre las líneas 400 y 490. 

Las sentencias IF-THEN que aparecen consiguen que no se pinten más puntos 
cuando una órbita se haya completado en el caso de los planetas y en el caso del 
cometa que no se pinten los puntos que quedan fuera de la pantalla, ya que 
saldría error. 

A partir de la línea 1130 se definen las ecuaciones del movimiento de los 
planetas, Halley y la sonda Giotto. Naturalmente representan las mismas elipses 
que las de las órbitas. 

La correspondiente al Giotto (x5, y5), partiendo de la Tierra, posee un radio 
a5 variable, para representar el movimiento espiral real que describe. 

La ecuación que define el movimiento del cometa no está parametrizada por el 
contador I como las demás, sino con el contador T. 

Esto tiene por objeto conseguir modular la velocidad del cometa con la 
variable Q, que tomará diferentes valores en cada tramo de órbita, para obtener 
la velocidad correcta del Halley en cada punto de su órbita. 

Además, el factor v0O de la línea 304 modula la velocidad de conjunto de todos 
los cuerpos; cada uno de estos tiene a su vez modulada su propioa velocidad por 
los argumentos de las sentencias SIN y COS. 

A diferencia de lo que sucedía con la versión MSX, aquí cada modelo de 
cometa está definido con un conjunto de puntos (plot). Cada uno de estos grupos 
entra en el momento adecuado en pantalla mediante las sentencias IF-THEN. 

Por ejemplo en las líneas 1520 y 1540 vemos que cuando XO esté entre 122 y 
105, y estemos en la mitad inferior de la pantalla (Y0*88) entrará en pantalla el 
cometa 4, con su velocidad correspondiente, y aparecerán sus características 
principales. 

Al igual que en el MSX, la sensación de movimiento se consigue con los 
contadores T e l. 

El borrado de cada conjunto de puntos (planeta, cometa o sonda) se consigue 
ejecutando dos veces las instrucciones (for n=0 to 1), la segunda de las cuales la 
lee con la sentencia OVER 1, invirtiendo el color de los puntos de la pantalla. 

La consecuencia de esto es que la segunda vez que lee la instrucción borra lo 
que había. El final del lazo n lo vemos en las líneas 1760, 1770 y 1780. El PAUSE 
100 consigue que estén más tiempo dibujados los cuerpos que borrados. 
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10 RETIRA A 


11 REM 

12 REM HALLE Y 

13 REM 

14 REM 44d IA 
15 REM 


30 PAPER Oz BORDER Oz INK 7: € 


31 PRINT AT 4,23" 


04d 


o TEOM LE 


CIRCLE 50,35,16 

37 PLOT 50,51: DRAW 180,-16 

38 PLOT 50,19: DRAW 180,16 

40 CIRCLE 50,35,14 

41 PLOT 50,45: DRAW 180,-10 

42 PLOT 50,25: DRAW 180,10 

50 CIRCLE 50,335, 12 

51 PLOT 50,39: DRAW 180,-4 

52 PLOT 50,31: DRAW 180,4 

60 CIRCLE 50,35,10 

70 CIRCLE 50,35,8 

80 CIRCLE 350,35,6 

90 CIRCLE 50,335,4 

100 CIRCLE 50,35,2 

280 PAUSE 200 

285 CLS 

290 PRINT "Este programa ha sid 
o pensado para representar lo 
mas fielmen—-te posible el encuen 
tro del co- meta Halley con la T 
ierra en 19896." 

292 PRINT "Todos los parametros 
son reales y las orbitas estan 
realizadas en perfecta escala." 
293 PRINT "La duracion del prog 
rama abarca doce meses.En cada t 
ramo apare- cen en pantalla las 
principales caracteristicas del 

cometa en ese momento.” 

294 PRINT "En el transcurso del 
programa se recoge la expedic 
ion de la sonda Giotto para el 
estudio delcometa Halley." 

300 PRINT : PRINT : PRINT 

301 PRINT TAB 9;3”"PARA COMENZAR” 
¿TAB 53"PULSE CUALQUIER TECLA” 


Ma, 
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302 IF INKEY*$="" THEN GO TO 30 


304 CLS : LET vo=2 

305 CIRCLE 115,88,1 

310 CIRCLE 115,89,2 

313. CIRCLE 113,889, 3 

317 BRIGHT O 

320 REM AI 

330 REM radios orbitas 

331 REM FITCH III EI ICE 

340 LET ai=25: REM mercurio 
350 LET a2=47: REM venus 

360 LET a3=66: REM tierra 

370 LET a4=100: REM marte 

380 LET a0=1180: LET bO0=296: RE 
M halley 

385 LET a5=aj-2: REM giotto 
390 REM orbitas planetarias 
392 LET e=.5 

393 LET mi=0: LET m2=0: LET n2= 
O: LET m3=0: LET m4=0 

394 LET ni=0: LET n3=0: LET n4= 
(8) 

393 FOR j=0 TO 9*PI STEP .17 
396 IF mi>143 AND ni1>86 AND j> 
i THEN GO TO 415 

397 REM AI III 

398 REM DIBUJAR ORBITAS 

399 REM III III EIA 

400 LET mi=120+a1+*C0S j 

410 LET ni=88+al*ex*SIN j 

415 IF m2>160 AND n2>88 AND j>1 
THEN GO TO 435 

420 LET m2=115+a2*C0S (3/2) 

430 LET n2=88+a2*ex*SIN (3/2) 
4353 IF m3>180 AND n3>87 AND j>2 
THEN GO TO 455 

440 LET m3=116+a3x*C0S (3/3) 

450 LET n3=88+a3*ex*SIN (3/3) 
4559 IF m4>222 AND n4>87 AND j>2 
THEN GO TO 480 

460 LET m4=124+a4x*C0OS (3/4) 

470 LET n4=88+a4*ex*SIN (3/4) 
480 LET mO0=1257+a0*C0OS (3/6) 
4790 LET nO0=88+b0*ex*SIN (3/6) 
1000 PLOT mi1,n1 
1010 PLOT m2,n2 
1020 PLOT m3,n3 
1030 PLOT m4,n4 
1034 IF m0>255 THEN GO TO 1050 
1033 IF n0>175 THEN GO TO 1050 
1036 IF noO THEN GO TO 1050 


1040 PLOT m0,n0 

1050 NEXT j 

1052 FOR m=0 TO 1 

10355 OVER 1: PRINT AT 0,03"2-7-8 
3": PRINT AT 19,0;3"El Giotto": P 
RINT AT 21,0;"sale al encuentro 
del Halley" : 
1060 IF m=1 THEN GO TO 1070 
1065 PAUSE 200 

1070 OVER Oz NEXT m 

1077 REM verde 0H IE MITRA 
1095 REM contadores; velocidades 
1097 REM 4d IE LISO IEA 
1100 LET I=-5 

1105 LET Q=0 

1107 LET T=0 

1110 LET I=1+5 

1113 LET a5=a5-.2 

1122 LET T=T+*Q0 

1123 REM 46-4-234-1636 3636-1636 26H 36-36 26-36 -4E363E 
1125 REM movimiento de planetas 
1127 REM 623636616 I0 IE II LE MITE MEL IEA 
1130 LET x1=120+a1*C0S (-2.92+.2 
1*1*v0/8.797) 

1140 LET yi=-88-exalx*SIN (-2.92+ 
.«21*1*v0/8.797) 

1150 LET x2=115+a2*C0S (-1.4+.21 
*I*v0/22.47) 

1150 LET y2=-889-ex*ra2x*SIN (-1.4+. 
Z21*1v0/22.47) 

1170 LET x3=116+a3*C0S (-2.8+.21 
*1*v0/36.526) 

1180 LET y3=-B88-e*tal3*SIN (-2.8+. 
21*1*v0/36.526) 

1183 LET x5=116+a5*C05 (-2.8+.21 
*1I*v0/26) 

11895 LET y5=88+a5*ex*SIN (-2.8+.2 
1*I*v0/26) 

1190 LET x4=124+a4x*CO0S (.21x*1*v0 
/69.6978) 

1200 LET y4=-B88-e*ra4*SIN (.21+*I1% 
vO/68. 679) 

1210 IF O=0 THEN  LET Q=5 

1229 REM ms ME IE IIA LE 

1230 REM HALLEY, ecuaciones 

1231 REM III IE MEE ETE RIE 

1240 LET x0=1257+2a0*C05 (2.5+T*v 
0/1000) 

1250 LET yO0=88-ex*rbO*SIN (2.5+TX*v 
0/1000) 

1260 IF yO0>165 THEN GOD TO 1740 
1270 IF y0>173 OR yO0<0 THEN LE 


YT y0=0 

1320 FOR n=0 TO 1: OVER 1 

1329 REM REE 

1330 REM planetas 

1331 REM dre 

1340 PLOT x1,y1 

1350 PLOT x1+1,y1 

1360 PLOT x1i+i,y1+1 

1370 PLOT x1,y1+1 

1420 PLOT x2,y2 

1421 PLOT x2,y2+1 

1424 PLOT x2,y2-1 

1426 PLOT x2+1,y2 

1427 PLOT x2+1,y2+1 

1430 PLOT x2+1,y2-1 

1438 PLOT x2-1,y2 

1439 PLOT x2-1,y2+1 

1441 PLOT x2-1,y2-1 

1442 PLOT x3,y3 

1443 PLOT x3,y3+1 

1444 PLOT x3,y3-1 

1445 PLOT x3+1,y3 

1446 PLOT x3+1,y3+1 

1447 PLOT x3+1,y3-1 

1448 PLOT x3-1,y3 

1449 PLOT x3-1,y3+1 

1450 PLOT x3-1.,y3-1 

1451 PLOT x4,y4 

1452 PLOT x4,y4+1 

1453 PLOT x4+1,y4 

1454 PLOT x4+1,y4+1 

1455 IF x5<x0 AND ySxXyO0 THEN GO 
TO 1470 

1458 PLOT x5,yS:z PLOT x5,yS-1:2 P 

LOT x35, y5-2 

1469 REM 4164-31-61 IE HARIAN AR 

1470 REM Halley, movimiento 

A 

1475 IF x0>255 THEN GO TO 1760 
1477 IF y0>88 THEN GO TO 1619 
1480 IF x0<200 AND y0<88 THEN 6 

O TO 1500 

1481 REM 4364232 

1482 REM cometa 1 

1483 REM A 

1485 LET 0=3.5 

1486 OVER Oz PRINT AT 0,03 "15-8- 

85/30-9-86": PRINT AT 2,03 "Magni 
tud=13/11": PRINT AT 18,03"D=19g 
.6m/20g.": PRINT AT 19,03”A.R.ó6h 
1m/6h12m": PRINT AT 21,0;3"Distan 
cia=550/315 millones de km": OVE 
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R 1 
1490 
1491 
1492 
1493 
1497 
1500 
501 
"1502 
1503 
1505 
1507 


85/24-10-85": 
itud= 
0./21g.12m": 


PLOT 
PLOT 


XO, yO 
xO+1,y0-1 
PLOT x0+2,yO-1 
PLOT x0+*3, y0-=1 
GO TO 17650 

IF x0<159 THEN 
REM 4636363424636 
REM cometa 2 
REM 43d 
LET Q=4 
OVER 0: 


GO TO 1520 


PRINT AT 0,03 "30-9-— 
PRINT AT 2,03 "Magn 
PRINT AT 18,03"D=2 

PRINT AT 19,03"A.R. 


11/9.5": 


£6hi2m/5Sh46m": PRINT AT 21,03 "Dis 
tancia=315/200 millones de km": 


OVER 
1510 
1511 
1512 
1513 
1514 
1517 
1520 
O TO 
1521 
1522 
1523 
1525 
1526 


—-83/15-12-86": 

nitud=9.5/6.2": 
-.R.5h46m/23h16m": 
"D=21g.12m/3g.37m": 


1 

PLOT 
PLOT 
PLOT 


x0, yO 
xO0+1,y0-—1 á 
x0+2, yO-1 a 
PLOT x0+2,y0-2 

PLOT 0+3, yD-2 

GO TO 1750 

IF x0<122 AND yO0<88 THEN 6G 
1540 

REM eee IR 

REM cometa 3 

REM 423144 

LET Q=5 ; 

OVER Oz PRINT AT 0,03"24-10 
PRINT AT 2,03 "Mag 
PRINT AT 19,0;3"A 
PRINT AT 18,05 
PRINT AT 21, 


0;3"Distancia=200/120 millones de 


km": 


1530 
1531 
1532 
1533 
1534 
1535 
1536 
1537 
1538 
1539 
1540 
O TO 
1541 
1542 
1543 
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OVER 1 

PLOT x0,y0 
PLOT x0+1,y0O-1 
PLOT x0+1,y0-2 
PLOT x0+3, y0-2 
PLOT x0+1,y0-3 
PLOT x0+3,y0-3 
PLOT x0+2,y0-4 
PLOT x0+4,y0-4 
PLOT x0+3,y0-5 
GO TO 1740 

IF x0<105 AND y0<88 THEN 6G 
1560 

REM IEA 
REM cometa 4 
o E 


Y 


1545 LET 0=6.5 

1546 OVER O: PRINT AT 0,03"15-12 
-85/07-01-86": PRINT AT 2,03"Mag 

nitud=6.2/5.6": PRINT AT 18,03"D 
=3g.37m/-3g.53m ": PRINT AT 19,0 
5"A.R.23h16m/22h2m PRINT AT 2 
1,03"Distancia=120/195 millones 
de km"z OVER 1 

1547 PLOT x0,y0 

1548 PLOT x0,y0-1 

15350 PLOT x0, y0-2 

13391 PLOT x0-1,y0-3 

1932 PLOT 021,03 

15593 PLOT x0,y0-4 

1554 PLOT x0-1,y0-5 

1555 PLOT x0+1,y0-5 

1559 GO TO 1760 

1560 IF x0<90 AND yO0<88 THEN GO 
TO 1580 
1561 REM 
1562 REM 


HERRIA 

cometa 3 

1563 REM Ad: 

1565 LET 0=9.5 

1566 OVER Oz PRINT AT 0,0;”7-1-8 
6/21-1-86 PRINT AT 2,05"Mag 
nitud=5.6/4.8": PRINT AT 18,0;"D 
=-3g.53m./-69.17m.  ": PRINT AT 
19,03"A4.R.22h2m/21h36m ": PRINT 
AT 21,03"Distancia=195/225 millo 
nes de km": OVER 1 

1568 PLOT x0,y0O 

1569 PLOT xO0-1,y0-1 

1570 PLOT x0-2,y0-2 

1371 PLOT 10-1,y0-3 

1372 PLOT x0-3,y0-3 

1573 PLOT x0-4,y0-4 

1574 PLOT x0-—2,y0-4 

1573 PLOT x0-3,y0-3 

1576 PLOT x0-4,y0-6 

1577 PLOT x0-2,y0-6 

1579 GO TO 1760 

1580 IF x0<79 THEN GO TO 1600 
1581 REM 443233336 

1582 REM cometa 64 

1583 REM A%41+% 343 

1585 LET Q=13.5 

1586 OVER Oz PRINT AT 0,03"21-O1 
-86/04-02-86": PRINT AT 2,0; "Mag 
nitud=4.8/4 ": PRINT AT 18,03"D 
=—69 . 28m/-—9g. 15m PRINT AT 19 
2053"A4.R.21h36m/2ih1im": PRINT AT 
21,03 "Distancia=225/230 millone 


s de km”: OVER 1 

1588 PLOT x0,y0 

1589 PLOT x0-1,y0 

1590 PLOT x0-4,y0-4 

1592 PLOT x0-1,y0-1 

1593 PLOT x0-3,y0-1 

1594 PLOT x0-2, y0-2 

1595 PLOT x0-4,y0-2 

1596 PLOT x0-3,y0-3 

1597 PLOT x0-6,y0-2 

1598 PLOT x0-5, yO0-3 

1599 GO TO 1740 

1500 IF x0>85 AND y0>87 THEN GO 
TO 1619 

1501 REM AHHH: 

1502 REM cometa 7 

1603 REM 40d 

1505 LET 0=17.5 

1606 OVER O: PRINT AT 0,035"04-02 

-86/15-02-86": PRINT AT 2,0; "Mag 

nitud=4/4.1 ": PRINT AT 18,05"D 

=-9g.13m/-11g.59%m ": PRINT AT 19 
105"A.R.21ih1im/20hS1m": PRINT AT 
21,03 "Distancia=230/2253 millone 
s de km": OVER 1 

1607 PLOT x0,y0 

1608 PLOT x0-2,y0 

16509 PLOT x0-3,yO0+1 

16510 PLOT x0-3,yO-1 

1611 PLOT x0-4,y0O 

1612 PLOT x0-5,y0+1 

1613 PLOT x0-S5, yO0-1 

1615 PLOT x0-6,y0 

1616 PLOT x0-7,y0-1 

1617 PLOT x0-7,y0+1 

1518 60 TO 1760 

1619 IF x0>92 AND y0>88 THEN GO 
TO 1440 

1520 LET 0Q=13.5 

16421 REM AMI 

1622 REM cometa 8 

FRANELA 

1623 IF x0>83S AND x0<88 THEN 
60 TO 1524 

1624 OVER O: PRINT AT 0,03"15-02 
—-86/01-03-86": PRINT AT 2,03"Mag 
nitud=4.1/4,4": PRINT AT 18,0;3"D 
=-11g.539m/—16g.19m":z PRINT AT 19 
,05"4.R.20h51m/20h27m": PRINT AT 
21,05 "Distancia=225/195 millone 
s de km": OVER 1 

1625 IF +0<85 OR x0>87 THEN GO 


TO 1627 
1626 OVER O: PRINT AT. 0,03"El Gi 
otto llega "z PRINT AT 2,0;3"a 
1 cometa ": PRINT AT 18,0; 
"...y se destruye “3 PRINT 
AT 19,03"al chocar con el ": PR 
INT AT 21,05"Sucede el 13 de Mar 
zo de 1986 ": PAUSE 100: OVER 1 
1627 PLOT x0,y0 
16528 PLOT x0-1,y0 
1629 PLOT x0-1,y0+1 
1630 PLOT x0-3,y0+1 
1631 PLOT x0-2, yO0+2 
1532 PLOT x0-4,y0+2 
1634 PLOT x0-3,yO0+3 
1635 PLOT x0-6, y0+2 
1636 PLOT x0-5, yO+3 
15637 PLOT x0-4,y0+4 
1639 GO TO 1740 
1540 IF x0>108 AND y0>88 THEN 6G 
O TO 1660 
1641 REM 64433 
1642 REM cometa 9 
1643 REM Are 
1544 LET 0=9.5 
1646 OVER Oz PRINT AT 0,03"01-03 
-86/21-03-86": PRINT AJ 2,03"Mag 
nitud=4.4/4.5": PRINT AT 18,03"D 
=-16g.19m/-26g.43m": PRINT AT 19 
.05"A.R.20h27m/19h38m”: PRINT AT 
21,0; "Distancia=195/120mi1llones 
de km": OVER 1 
16547 PLOT x0,y0O 
1548 PLOT x0-1,yO+1 
1549 PLOT x0-2, yO+2 
1550 PLOT x0-1,y0+3 
1651 PLOT x0-3,yO+3 
1552 PLOT x0-4, yO0+4 
15653 PLOT 0-2, yO+4 
1654 PLOT x0-3,y0+5 
1655 PLOT x0-4,yO+6 
1556 PLOT x0-2, y0+6 
1559 GO TO 1740 
1560 IF x0>124 AND y0>88 THEN G 
O TO 14580 
1661 REM 436442631 
1662 REM cometa 10 
1663 REM MH -46 2H AE 
1665 LET 0O=6.5 
1666 OVER O: PRINT AT 0,03"21-03 
-86/09-04-86": PRINT AT 2,0; "Mag 
nitud=4.5/3.9": PRINT AT 18,0;3"D 
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=-269.43m/-479q.23m":2 PRINT AT 19 
103"A.R.19h38m/15h48m": PRINT AT 
21,05"Distancia=120/70millones 

de km": OVER 1 

1567 PLOT x0,y0 

1568 PLOT x0,yO+1 

1669 PLOT x0-1,yO0+2 

1670 PLOT xO+1,yO0+2 

1571 PLOT x0,y0+3 

1572 PLOT x0-—1,y0+4 

1573 PLOT x0+1,y0+4 

1574 PLOT x0,yO0+5 

16575 PLOT x0-—1,y0+6 

1576 PLOT x0+1,y0+6 

1579 GO TO 1740 

1580 IF x0>159 AND yO0>98 THEN G 
Q TO 1700 

156891 REM 3334-4234 

1582 REM cometa 11 

1583 REM 4333 

15835 LET 0=5 

1586 OVER O: PRINT AT 0,03"09-04 
—86/08-05-86": PRINT AT 2,03 "Mag 
nitud=3.9/6.8": PRINT AT 18,03"D 
=-473.23m/-13qg.26m": PRINT AT 19 
.03"A.R.15h22m/10h38m": PRINT AT 
21,0; "Distancia=70/150 millones 
de km "z OVER 1 

1587 PLOT x0,y0 

1688 PLOT x0,yO+1 

15899 PLOT x0,y0+2 

1590 PLOT x0+1,y0+1 

1591 PLOT x0+2,y0+2 

1592 PLOT x0+1,y0+3 

1693 PLOT x0+3,y0O+3 

1594 PLOT x0+2,y0+4 

1695 PLOT x0O+1,y0+5 

1599 GO TO 17650 

1700 IF x<0>200 AND yO>88 THEN 6G 

O TO 1720 

1701 REM eee 

1702 REM cometa 12 

1703 REM 3463332363636 26 

17035 LET 0=4 

1706 OVER Oz PRINT AT 0,03"08-05 
—86/15-06-86": PRINT AT 2,03 "Mag 
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nitud=6.9/9.4": PRINT AT 18,05"D 
=—-1i3g.26m/-59.21m ":z PRINT AT*19 
,05"A.R.10h38m/10h26m": PRINT AT 
21,03 "Distancia=150/330 millone 
s de km": OVER 1 

1707 PLOT x0,y0O 

1708 PLOT x0+1,y0+1 

1709 PLOT x0+2,y0+1 

1710. PLOT x0+2, yO+2 

1711 PLOT x0+3, yO0+2 

1719 GO TOD 1760 

1721 REM AI 

1722 REM cometa 13 

1723 REM 1433533336363 

1725 LET Q=3.5 

1726 OVER O: PRINT AT 0,05" 15-06 

—86/15-08-86": PRINT AT 2,05”Mag 

nitud=9.4/11.6": PRINT AT 18,03” 
D=-5g.21im/-69.43m  "z PRINT AT 1 
9,05"AM.R.10h26m/11h0m ": PRINT A 

Y 21,03”"Distancia=330/570 millon 
es de km": OVER 1 

1727 PLOT x0,y0 

1728 PLOT x0+1,y0+1 

1729 PLOT x0+2,yO+1 

1730 PLOT x0+3,yO+1 

1740 IF x0>255 THEN GOD TO 5000 
1760 IF n=1 THEN GO TO 1780 
1770 PAUSE 100 

1780 OVER Oz NEXT n 

1800 GO TO 1110 

5000 CLS 

5010 PRINT "Este ha sido el reco 

rrido del cometa Halley en su 
acercamien- to a la Tierra en 19 

856." 

5015 PRINT : PRINT : PRINT : PRI 

NT 

5020 PRINT TAB 2;3"Para volver a 

empezar pulsar"3TAB 853 "cual qui er 
tecla." 

5023 PAUSE 100 

5030 IF INKEY*$="" THEN GO TO 50 

30 

5040 CLS : GO TO 30 


sr. 


Hagnrnitud=4,., 3L 


LD DEA 


ens Io 


Hagnitud=z:d. 1-4. ¿4 


LLanes 


Figura 10. Volcados de pantalla de la versión para Spectrum del programa “Halley” 
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Programa número b: El Halley sobre el cielo 
nocturno (versión MSX) 


Este programa permite localizar al Halley en el fondo de la bóveda celeste. Se 
representa la porción de cielo que “recorrerá” aparentemente. 

Las relaciones de tamaño y forma de las constelaciones son las reales, ya que 
se han dibujado a escala (véase Fig. 11). 


30 PLEIADES. A PEGASUS 
z Nov 1 
2 20 
Q . Nov21 AQUILA  OPHIUCUS 
Zz 
e ici REGULUS 
a.0 
30 
-20 
CAPRICORNUS ' ANTARES 
pa . e, HYDRA 
FOMALHAUT dai 


CENTAURUS 


6 4 2 24 22 20 18 16 14 12 10 
ASCENSION RECTA 


Figura 11. Constelaciones del hemisferio norte y recorrido del cometa “Halley” 


Los diferentes sprites del Halley representan la morfología con la que lo 
observarmos, pero exagerando algo el tamaño relativo con el fondo estrellado 
para su mejor apreciación. 

Antes del perihelio, la cola se extiende hacia el Este. A medida que la cola gire 
respecto a la Tierra, veremos al cometa más “de frente”, es decir, la cola irá 
aparentemente disminuyendo hasta “desaparecer” (junto con el cometa) en el 
perihelio, al situarse “detrás del Sol”, aunque en la realidad sea el momento de 
mayor actividad. Pasado el perihelio, ocurre lo contrario: la cola aumenta de 
tamaño aparente (aunque en la realidad disminuye). 

Las diversas posiciones del cometa sobre el fondo y las previsiones en cuanto 
a su morfología son las previstas, si no coinciden exactamente con la realidad 
pueden modificarse los valores en el listado. 
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Detalles técnicos del programa 


Desde el principio del listado hasta la línea 730 se dibujan las constelaciones 
sobre un fondo oscuro. Las estrellas se representan mediante puntos (PSET), 
cuyas coordenadas son leídas mediante un bucle por las sentencias READ-DATA. 

Mediante dos sentencias CIRCLE en las líneas 990 y 1000, se representa la 
trayectoria del Halley, aproximada a su recorrido y válida para la resolución de la 
pantalla. 

Para trazar los números cercanos a la trayectoria se utilizan sentencias LINE 
a partir de la línea 1020. Las características del Halley al pasar por esos seis 
puntos aparecen en la parte inferior de la pantalla. Dado que estamos utilizando 
una pantalla de gráficos, es necesario abrir un fichero para escribir estos datos, 
línea 1170. 

Las ecuaciones del movimiento del Halley se definen en las líneas 1380 y 1550. 
Como hemos aproximado la trayectoria mediante dos curvas, necesitaremos dos 
pares de coordenadas (X1, Y1) y (X2, Y2), parametrizadas por las variables L y K. 

R y S con las variables que determinan el sprite que debe de utilizarse en cada 
momento, mediante las sentencias IF-THEN de las líneas 1390 a 1550. 

El RESTORE de la línea 1660 nos permite volver a leer los DATA, si 
queremos volver a ejecutar el programa. 

La velocidad del cometa está determinada en la línea 1350 por la variable V. 


A Y 

2 * 160 ” LAZO DE LECTURA DE DATOS 

3 ?”% MOVIMIENTO DEL HALLEY SOBRE + A 
* EL CIELO NOCTURNO * 180 > 

e E + 190 READ A,B 

RR III MEM EMINEM LIM IMHE 200 PSET (A+20,B+5),15 

6? 210 NEXT I 

10 KEY OFF:COLOR 15,1,1 220 > 

20 CLS 230” AAA IM MAITE MDI HET EE IDR ME HERE 

30 CLEAR 


40 SCREEN O 

50 PRINT"Este programa intenta dar un 
a vision aproximada del aspecto del c 
ometa Halley sobre el cielo noctur 
no en su aproximacion a la Tierra de 

1984." 

60 PRINT 

70 PRINT"Comprende desde el 25 de Nov 
iembre del985 hasta el 30 de Mayo de 

1996, coincide,por tanto,con el pe 
riodo en el que es observable a simpl 
e vista, excepto a su paso por el per 
ihelio.” 

80 LOCATE 0,20 

90 PRINT"Para comenzar pulse barra es 
paciadora"” 

100 IF STRIG(0)=0 THEN 100 

110 CLS 

120 SCREEN 2,2 

130 FOR 1=0 TO 104 

140 ” 


240 ? POSICIONES DE LAS ESTRELLAS 
DE CADA CONSTELACION 

250 ARA 

260 ” 

270 ” TRIANGULO 

280 >” 

290 DATA 14,4,10,7,18,9 


310 ”>PEGASO 


330 DATA 38,10,36,21,54,10,69,12,61,1 
4,53,23,62,29,67,32,72,30 

340 ” 

350 *CETUS (BALLENA) 

360 ” 

370 DATA 2,32,6,45,10,27,17,32,17,42, 
18, 48,26,54,32,62,42,57,45,58,48,65,3 
5,66,32,74,1,73,15,94 

380 ” 

390 ” CAPRICORNIO 

400 ” 

410 DATA 73,64,75,657,84,64,93,59,92,6 
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1,81,78,90,82,69,87 

420 ” 

430 ” AGUILA 

440 > 

450 DATA 99,21,96,23,89,36,99,35,105, 


460 ” 

470 ”>SAGITARIO 

480 ” 

490 DATA 111,77,110,72,115,72,117,73, 
120,74,122,67,122,82,124,898,127,83 
500 ” 

510 ”>OPHIUCHUS 

5z0 ? 

530 DATA 127,22,129,14,134,17,142,23, 
143,27,141,28,140,33,131,36 

540 ” 

550 ”ESCORPION 

560 ” 

570 DATA 141,92,138,91,136,95,138,98, 
141,99,145,97,148,98,149,92,151,82,15 
5,78,161,77,159,679,156,67,163,85,163, 
88 

580 ?” 

590 ”VIRGO 

£00 ” 

610 DATA 158,45, 163, 44,174,47,177,42, 
172,36,167,32,183,38,187,36,182,32,18 
3,24,178,38 

620 ” 

630 ” CENTAURO 

640 ”? 

650 DATA 184,94,191,97,196,105,193,12 
2,202,97,198,86,189,87 

660 ? 

670 “HYDRA 

680 ?” 

690 DATA 196,71,205,66,217,72,220,71, 
223,56,225, 58,229,534 

700 ” 

710 *LEO 

720 ” 

730 DATA 202,22,207,13,210,24,220,29, 
223,25,222,179,225,16,205,32 

740 ” 

7530 ” 

A 
770 ” DEFINICION DE LOS SPRITES 
O 
790 ” 

800 'A1$=CHR$ (LYHO) +CHR$ (HO) +CHR$ (LHO) 
+CHR$ (HO) +CHR$ (HO) +CHR$ (HO) +FCHRS$ (Ze 
HO) +CHR$ (HO) 

810 C1$=CHR$ (HO) +CHR$ (LHO ) FEHR$ (LHO) 
+CHR$ (HO) +CHR$ (2H8) +CHR$ (2H 1 4) +CHRS ( 
HA) +CHRS$ (%:H1 4) 

820 SPRITES (1)=C15$ 

830 C2%$=CHR$ (HO) +CHR$ (HO) +CHR$ (HO) 
+CHR$ (HO) +CHR$ (%H1 4) +¿CHR$ (2%2H28) +CHR$ 
(2H15)+CHR$ (2H3) 

840 SPRITES (2) =C2$ 
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850 C3$=CHR$ (KHO) +CHR$ (RHO ) +FCHRS$ (%H10O 
) +CHR$ (2428) +CHR$ (%H54 ) +CHR$ (2H28) +CH 
R$ (%H16) +CHR$ (8H3) 

860 SPRITES (3) =C3% 

870 C4$=CHR$ (LHCO) +CHR$ (%H68) +CHR$ (3H 
14) +CHR$ (2H28) +CHR$ (%HO) +CHR$ (HO) +CH 
R$ (HO) +CHR$ (LHO) 

880 SPRITES (4) =C4$ 

890 CS$=+CHR$ (LHCO) +CHR$ (2H568) +CHR$ (2% 
H14)+CHR$ (%H24) +CHRS (4H14) +CHR$ (Y H89) + 
CHR$ (2:H0) +CHR$ (2%HO) 

900 SPRITES (5) =C5$ 

910 Có$=CHR$ (2H28) +CHR$ (%HSO) +CHR$ (2H 
28) +CHR$ (4H10) +CHR$ (HO) +CHR$ (HO) +CH 
R$ (HO) +FCHR$ (HO) 

920 SPRITES (6) =C6$ 

930 SPRITES (0)=A1%+ 


E AAA 
960 ” TRAYECTORIA Y REFERENCIAS 
A AAA 
980 ” 

990 CIRCLE (-30,215),203,4,.7,1. 
1000 CIRCLE (177,45),65,4,3.9,6. 
1010 COLOR11 

1020 LINE (20,27)-(20,22) 

1030 LINE (70,45)-(70,50) 

1040 LINE (67,45)- (67,50) 

1050 LINE (97,62)- 197,67) 

1060 LINE (94,62)- (94,67) 

1070 LINE (91,62)- (91,67) 

1080 LINE (135, 105)-(135,110) 
1090 LINE (138, 105)- (140,110) 
1100 LINE (140,110)-(142,105) 
1110 LINE (180,115)-(182,120) 
1120 LINE (182,120)-(194,115) 
1130 LINE (246,45)- (246,40) 

1140 LINE (239,40)-(241,45) 

1150 LINE (243,40)-(241,45) 

1150 COLOR1S 

1170 OPEN"GRP: "FOR OUTPUT AS+1 
1180 PRESET (10, 129) 

1190 PRINT+1,"I :25-11-28553Mg ó6jEn A 
ries" 

1200 PRESET (10,140) 

1210 PRINTH41,"11 :15-12-853Mg -5.85En 
Piscis” 

1220 PRESET (10,151) 

1230 PRINT+1,"111:9-2-863perihelio;zNo 
visible” 

1240 PRESET (10, 162) 

1250 PRINTH1,"IV :25-3-2863Mg 45En Sag 
itario" 

1260 PRESET (10,173) 

1270 PRINT+1,"W  :10-4-863Mg 4.23En C 

entauro" 

1280 PRESET (10, 184) 

12790 PRINT+1i,”"VI :30-5-863Mg 653En Hyd 

ra" 

1300 ” 

A E 


6 
17 


úl CONTADOR Y POSICIONES 
PANA RIIIE INMI IAEA 
V=1 
=14+4.9 
T=T+5E-O3*V 
X1=-30+203*C0S (L.) : Y1=2135+203*SIN 


IF X1>133 THEN 1480 


IF X1i<-10 OR X1>103 THEN R=0 
IF X1<40 AND X1>-10 THEN R=1 
IF X1<60AND X1>40 THEN R=2 

IF Xi<80 AND X1>60 THEN R=3 
IF X1<103 AND X1>80 THEN R=2 
IF X1>115 THEN R=4 

PUT SPRITE 4, (X1-8,Y1-8),15,R 
GOTO 1360 

K=I1+35.63 

I=1+.015*VY 
X2=177+65*C0S (K) : Y2=45-65*S IN (K) 
IF X2<145 THEN S=4 


1520 IF X2>145AND X2<180 THEN S=5 
1530 IF X2>1804ND X2<215 THEN S=4 
1540 IF X2>215AND X2<240. THEN S=6 
1550 IF X2>240 THEN S=0 

1560 PUT SPRITE 3, (X2,Y2),,S 

1570 PUT SPRITE 4, (X1-9,Y1-9),,0 

1580 IF X2>241THEN 1600 

1590 GOTO 1480 

1600 CLS 

1610 SCREEN O 

1520 PRINT"Este ha sido el recorrido 
aproximado del cometa Halley hasta +fi 
nales de Mayo de 1986,momento a par 
tir del cual deja de ser posible s 
u observa- cion a simple vista." 

1630 LOCATE 1,20 

1540 PRINT"Si quiere volver a empezar 
pulse la barra espaciadora." 

1450 IF STRIG(0)=0 THEN 1650 

1660 RESTORE 

16570 GOTO 20 
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El 


Halley sobre el cielo nocturno (versión Spectrum) 


Este programa intenta representar la trayectoria que seguirá el Halley sobre el 
cielo nocturno, durante el período en que será visible. 
En él se incluyen los datos más importantes del cometa en algunos puntos de 


su órbita. 


Primero hacemos los círculos que representarán la órbita del Halley (líneas 
250-320). Mediante los PLOT volvemos a hacer los diferentes diseños del cometa 
que irán entrando a lo largo de la trayectoria. 

Con las sentencias IF-THEN de las líneas 1140 a 1200, hacemos entrar cada 
uno de ellos en el momento adecuado en el primer tramo de la órbita (antes del 


perihelio). 


Y con las de las líneas 1440 a 2120, lo mismo para la segunda parte de la 


órbita (después del perihelio). 


Las ecuaciones de movimiento del cometa están en las líneas 1120 para el 
primer tramo y 1400 para el segundo. 
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Figura 12. Volcado de pantalla del programa. “El Halley sobre el cielo nocturno” (versión 
Spectrum) 


REM Hem IE III II MEE 

REM EL HALLEY SOBRE EL 
CIELO NOCTURNO 

REM AIM EA II II MEM III 

cLs 

BORDER Oz PAPER O: INK 7 

FOR I=0 TO 104 


220 
230 
240 
241 
243 
245 
250 


READ A,B 

PLOT A+20,-B+175 

NEXT 1 

REM Id IEEE DEJE FEE DE ME HE IE EE 
REM DIBUJO DE LA ORBITA 
REM III AI MEME IET TER DEERE EE 
FOR j¡=0 TO .86 STEP .02 


1] 

2559 LET mi=-30-203*C0S (j+1.6): LET 
n1=-40+203*SIN (j+1.6) 

260 IF mi>=255 THEN GO TO 290 

263 IF mi<=0 THEN GO TO 290 

266 I1F n1i<=0 THEN GO TO 290 

269 1F ni>=173 THEN 60 TO 290 

275 PLOT m1,ni 

290 NEXT j 

295 FOR i=0 TO .8*PI STEP .06 

300 LET m2=177+46435*C0S (1+3.8): 

LET n2=130+65*SIN (i+3.8) 

305 IF m2<=0 THEN 60 TO 320 

308 IF m2>=255 THEN GO TO 320 

311 IF n2<=0 THEN (GO TO 320 

314 IF n2>=1735 THEN GO TO 320 

317 PLOT m2,n2 

320 NEXT i 

325 REM 3H DE IE IEEE 

327 REM DATOS DEL HALLEY 

330 REM 43M IE MIE MIER II 

340 PRINT AT 16,053"1I 25-11-853Mg 6jE 
n Aries" 

3530 PRINT AT 17,03"II 15-12-855Mg 5. 
8;En Piscis" 

360 PRINT AT 18,0;"111 9-2-B6;5Perihe 
lios5No visible" 

370 PRINT AT 19,03"IV 26-3-863Mg 45E 
n Sagitario” 

380 PRINT AT 20,0;”"V 10-4-B63Mg 4.2; 
En Centauro" 

390 PRINT AT 21,0;"VI 30-5-B6;Mg 6;3E 
n Hydra" 

890 PRINT AT 3,2;3”I" 

900 PRINT AT 4,83"11" 

910 PRINT AT 8,103"III" 

950 PRINT AT 13,163"IV" 

990 PRINT AT 14,233"vV" 

1000 PRINT AT 5,293 "VI" 

1105 REM 43360363636 HI IEA IRM RAR 

1106 REM ECUACIONES DEL HALLEY 

A E 

1110 FOR t=0 TO .87 STEP .03 ñ 

1120 LET x1=-30-203xCOS (t+1.5): LET 
y 1=-40+203*SIN (t+1.5) 

1130 LET x1i=x1+1: LET yi=y1+1 

1140 IF x1>133 THEN (GO TO 1380 

1150 IF x1<2 THEN GO TO 1350 

1165 FOR n=0 TO 1: OVER 1 

1168 IF x1>115 THEN G0 TO 1250 

1170 IF x1>103 THEN GO TO 1279 

1180 IF x1>80 THEN GO TO 1230 

1190 IF x1>60 THEN GO TO 1240 

1200 IF x1>40 THEN GO TO 1230 

1220 PLOT x1,y1z PLOT x1,y1+1: PLOT x 
1-1,y1: PLOT x1-1,y1+2: PLOT x1-2,y1+ 
1: GO TO 1280 

1230 PLOT x1,y1: PLOT x1-1,y1+1: PLOT 
x1-2,y1+22 PLOT x1-3,y1+32 PLOT x1-3 
syi+r1: PLOT x1-4,y1+2: GO TO 1280 
1240 PLOT x1,y1:z PLOT xi-1,y1:z PLOT x 


1i-1,yi+1: PLOT x1-2,y1+2: PLOT x1i-3,y 
1+3: PLOT x1-4,y1+4: PLOT x1-2,y1+1: 

PLOT x1-3,y1+1: PLOT x1-4,y1+2: PLOT 

Xx1-5,y1+3: PLOT x1-6,y1+2: PLOT x1-5, 

yi+1: GO TO 1280 

1250 PLOT x1,yl1:z PLOT xi+1,y1-1:3 PLOT 
x1+2,y1-2: PLOT x1+3,y1-1:2 PLOT x1+3 
sy1-3: PLOT x1+4,y1-2: GO TO 1280 
1279 PAUSE 25 

1280 IF n=1 THEN GO TO 1300 

1290 PAUSE 80 

1300 OVER Oz NEXT n 

1350 NEXT t 

1380 FOR k=3.8 TO 1000 STEP .1 

1400 LET x2=177+65*C0S k: LET y2=130+ 
65*SIN k 

1410 LET y2=y2+1  . 

1440 IF x2>240 THEN GO TO 2300 

2065 FOR m=0 TO 1: OVER 1 

2080 IF x2>215 THEN GO TO 2150 

2090 IF x2>180 THEN GO TO 2120 

2100 IF x2>145 THEN GO TO 2130 

2120 PLOT x2,y2: PLOT x2+1,y2-1: PLOT 
x2+2,y2-22 PLOT x2+3,y2-1: PLOT x2+3 
sy2-3: PLOT x2+4,y2-2: GO TO 2180 
2130 PLOT x2,y2: PLOT x2+1,y2: PLOT x 
2+1,y2-1: PLOT x2+2,y2-1: PLOT x2+2,y 
2-2: PLOT x2+3,y2-1: PLOT x2+3,y2-3: 
PLOT x2+4,y2-2:3 PLOT x2+4,y2-4: PLOT 
x2+5,y2-1: PLOT x2+5,y2-3: PLOT x2+6, 
y2-2: GD TO 2180 

2150 PLOT x2,y2: PLOT x2,y2+1: PLOT x 
2+1,y2: PLOT x2+1,y2+2:3 PLOT x2+2, y2+ 
1: GO TO 2180 

2180 IF m=1 THEN GO TO 2200 


2400 PRINT "SI QUIERE VOLVER A EMPEZA 
R,PULSE CUALQUIER TECLA” 

2450 IF INKEY$="" THEN GOD TO 2450 
2470 RESTORE 

25300 GO TO 10 

25360 REM 61M III IRA 

2545 REM CONSTELACIONES 

2570 REM III IE ME 

2600 REM Triangulo 

2700 DATA 14,4,10,7,19,9 

2800 REM Pegaso 

2900 DATA 38,10,36,21,54,10,69,12,61, 
14,53,23,62,28,067,32,72,30 

2950 REM Cetus 

3000 DATA 2,32,6,45,10,27,17,32,17,42 
,189,48,26,54,32,62,42,57,45,58, 48, 65, 
353,166,32,74,1,73,15,94 

3100 REM Capricornio 

3200 DATA 73,66,75,67,84,64,93,579,92, 
61,81,789,70,82,69,87 

3300 REM Aguila 


79 


3400 DATA 99,21,796,23,89,36,97,35,105 
. 59 

3500 REM Sagitario 

3600 DATA 111,77,110,72,115,72,117,75 
y120,74,122,67,122,82,124,88,127,83 
3700 REM Ofiuco 

3800 DATA 127,22,129,14,134,17,142,23 
,143,27,141,28,140,33,131,36 

3900 REM Escorpion 

4000 DATA 141,92,138,791,136,95,138,98 
141,99, 145,97,148,88, 149,92,151,82,1 
55,78,161,77,159,69,156,67,163,85,165 
, 88 


$ 
4100 REM Virgo 
4200 DATA 158,45,163,44,174,47,177,42 
,172,36,167,32,183,38,187,36,182,32,1 
83,24,1978,38 
4300 REM Centauro 
4400 DATA 184,94,191,97,1796,105,193,1 
22,202,97,198,86, 189,87 
4500 REM Hydra 
4600 DATA 196,71,205,66,217,72,220,71 
$ 223,56, 2259,59,229,54 
4700 REM Leo 
4800 DATA 202,22,207,13,210,24,220,29 
$223, 29,222+19,2239,16,205,32 


Apéndice A 


Realización de 
programas astronómicos 
en BASIC para MSX 


Los siguientes programas numéricos nos permitirán conocer rápidamente, en 
el primero, la magnitud aparente prevista para el Halley y, el segundo, la altura 
sobre el horizonte de cualquier astro según el lugar y la hora de observación. 

Será recomendable conocer estos datos si queremos observar el Halley, en 
especial cuando, debido a su poca luminosidad, sea necesario un pequeño 
instrumento para su localización. 

Estos son algunos procedimientos útiles para hacer programas astronómicos 
en BASIC de MSX. 

Para trazar órbitas se utilizará la sentencia CIRCLE con la excentricidad 
adecuada. 

Para representar astros se usarán sprites. Un buen método es dibujarlos 
previamente en un papel cuadriculado y buscar posteriormente el código corres- 
pondiente a cada línea. Recuérdese que definiendo cuatro bloques de 8 x 8 
puntos se puede construir un sprite de 16 x 16. A la hora de listar la definición de 
los sprites se ahorrará tiempo designando a una tecla de función la secuencia 
) + CHRS(8H”. 

Para describir trayectorias elípicas, la parametrización más sencilla es: 


coordenada X=A x cos(T) 
coordenada Y=B x sin(T) 


donde T es una variable que se incrementa lentamente entre O y 2 z (pi...) y A, B 
son los semiejes de la elipse. 
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Para mover un sprite utilizaremos una sentencia PUT SPRITE en el punto 
X, Y. A medida que aumenta el valor de T, éste nos describirá la trayectoria 
elíptica. 

Cuando queramos mover simultáneamente más de un sprite a diferentes 
velocidades, podremos utilizar varios contadores, o uno solo introduciendo 
factores que modifiquen el argumento de las funciones seno y coseno (en este 
caso T). 

Si queremos modular la velocidad del movimiento, haremos que el incremento 
de T sea variable, según la función deseada. Si, por ejemplo, hacemos que esa 
función sea inversamente proporcional al cuadrado de la distancia al foco de la 
elipse (como hicimos en los programas 1 y 3), la órbita resultante seguirá la 


segunda ley de Kepler. 


El punto de comienzo del movimiento se consigue con un desfase, dando un 


valor inicial constante a T. 


Para escribir caracteres en una pantalla de gráficos, será necesario abrir un 
fichero considerando la pantalla como un periférico, mediante la sentencia OPEN. 

Para aumentar el área de memoria de caracteres y poner a cero todas las 
variables, se utiliza la sentencia CLEAR. 

Si un programa va a ejecutarse repetidamente, no se debe olvidar utilizar la 
sentencia CLEAR, pues de lo contrario el valor inicial de las variables será el final 
de la anterior ejecución. Asimismo, se debe utilizar la sentencia RESTORE para 


volver a leer los DATA. 


Un fondo estrellado es fácil de conseguir con la sentencia PSET dentro de un 
bucle de lectura (READ) de las coordenadas en una sentencia DATA, según se ha 


hecho en el quinto programa. 


E 
2 >x E 
3% CALCULO APROXIMADO DE e 
4 ”* LA MAGNITUD TOTAL APARENTE + 
DD +. DEL COMETA HALLEY * 
6 "+ * 
7 


PHN NE MIMI MI MIRE 

8 ES 

10 KEYOFF : CLS 

20 PRINT "A partir de la distancia de 
l cometa ala Tierra y al Sol,expres 
ada en u- nidades astron"micas de dis 
tancia (u.a.31 u.a.= 150 millones 

de Km.), podemos calcular la magnitu 

d aproxi-— mada que tendra el Halley." 
30 PRINT 

40 PRINT "Los calculos se realizan de 
acuerdo con las formulas deducidas 

por Morrisy Green." 

50 PRINT:PRINT 

60 INPUT "Distancia aproximada a la T 
ierra (En unidades astronomicas)" 
5A 

70 IF A<1 THEN 60 

80 PRINT:PRINT 

90 INPUT "Distancia aproximada al Sol 

(En unidades astronomicas)" 
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100 IF B<1 THEN 90 

110 X=LOG(A)/LOG(10) : Y=LOG (B) /LOG (10) 

120 PRINT:PRINT 

130 INPUT "Antes del perihelio:A; 

despues de el:D"3A5$ 

140 IF A$="D" OR A$="d" THEN 190 

150 IF A$="A" OR Ab$="a" THEN 170 

160 GOTO 130 

170 M=5.47+5*X+11.1x*Y 

180 GOTO 200 

190 M=3. 14+5%X4+7. 75% Y 

200 PRINT:PRINT 

210 CLS 

220 LOCATE 0,8 

230 PRINT "desd 16d EI MIE E DIE REAEME IEA MIRA 

FRIA? 

240 PRINT USING ” 

MADA :HH.4H"5M 

250 PRINT 

260 PRINT "de ve-2-16-366 16-36-36 16H HEN HE IE MIMI 

MANERA 

270 LOCATE 2,20 

280 PRINT"Para volver a empezar pulse 
barra":PRINT TAB(13)5 "espaciadora" 

290 IF STRIG(0)=0 THEN 290 

300 RUN 


MAGNITUD APROXI 


7 HEM MEM RIES SIR MINIMA 
> + * 
? * CALCULO DE LA ALTURA + 
E DEL HALLEY * 
>» * 
O e 
> 


NO 0d 


10 KEYOFF: CLS 

20 LOCATE 0,3 

30 PRINT"Este programa es valido para 
calcularla altura sobre el horizonte 
en gra- dos de cualquier astro, aunqu 
e todos los comentarios esten referi 
dos al cometa Halley.” 

40 LOCATE 5,20 

50 PRINT "Pulse la barra espaciadora. 
60 IF STRIG(0)=0 THEN 40 

70 CELS 

80 PRINT "Uno de los datos necesarios 
para cal-—-cular la altura del Halley 
sobre el horizonte es el tiempo side 
reo local.La forma de calcularlo es l 
a siguien—-te:":PRINT 

90 PRINT "ix Contar el numero de mese 
s enteros transcurridos desde el 22 d 
e Septim- bre anterior y multiplicarl 

os por 2.” 

100 PRINT 

110 PRINT "2x Sumar la hora local (la 
oficial-1,o -2 en verano) del observa 
dor a este resultado.Esto es la hora 
siderea lo-cal (si fuera superior a 24 
restar 24)" 

120 PRINT "3x Los minutos se obtienen 
contando el numero de dias enteros 
transcurri-dos desde el 22 del mes an 
terior y multiplicandolos por 4 (lo 
gicamente, cada 60 minutos seran 1 ho 
ra mas a sumar al anterior resultad 

o.” 

130 PRINT 

140 PRINT TAB(7)3"Pulse barra espacia 
dora" 

150 IF STRIG(0)=0 THEN 150 

160 CLS 

170 INPUT "LATITUD DEL OBSERVADOR, gra 


dos "A 

180 INPUT "minutos "3B 

190 PRINT 

200 INPUT "LONGITUD DEL OBSERVADOR, gr 
ados"53D 

210 INPUT "minutos”;5 

220 PRINT 

230 INPUT "COORDENADAS DEL ASTRO, Decl 
inacion,grados "5; 

240 INPUT "minutos "3K 


260 INPUT "Ascension recta,horas "36 
270 INPUT "minutos "5 


290 INPUT "TIEMPO SIDEREO LOCAL,horas 
300 INPUT "minutos"3N 


320 INPUT "HORA DE OBSERVACION, horas” 


330 0=(B/560)+A 

340 S=(H/60)+6 

350 T=1(K/60)+3J 

360 U=(N/60)-+M 

370 R=(E/60)+D 

380 V=1.0027379094P+U+R 

390 X=V-S5 

400 Y=SIN(T)%SIN (0) +C0S (T) COS (Q) *COS 


410 Z=SQR (1-Yx+*Y) 

420 W=Y/Z 

430 CLS 

440 LOCATE 0,7 

A 

MAMI IE MAH 

460 PRINT 

470 PRINT USING " ALTURA DEL HALLEY ( 

en grados) H4t+. 4" 5 ATN (W) *180/3. 1416 

480 PRINT 

490 PRINT "04001066 HE IE IM MIMI 

RARA PRINT 

500 LOCATE 0,20 

510 PRINT"Para volver a empezar pulse 
la barra":PRINT TAB(13)3 "espaciadora 
s 


520 IF STRIG(0)=0 THEN 520 
330 GOTO 10 


Apéndice B 


Cálculo de la posición 
y parámetros del Halley 


REM 
REM 
REM 
REM 
REM 
REM 
REM 


en el IBM PC 


Se añade a continuación el listado de un programa escrito en BASICA sobre 
un IBM PC que calcula los parámetros básicos de la trayectoria aparente del 
cometa con la intención de que pueda ser utilizado por los usuarios de compati- 
bles PC. Por otra parte, el programa es fácilmente traducible a cualquier BASIC 
que tenga la precisión suficiente (más de diez dígitos por número). 

Hay que hacer notar que los datos generados por el programa sólo son 
válidos cuando el cometa se halla cerca del perihelio de la órbita, esto es, tres 
meses por delante y por detrás del momento en que se produce éste, el 9 de 
febrero de 1986. 

El programa es una aplicación general; por tanto, se puede aplicar al estudio 
de cualquier cometa del que se conozcan los datos de su órbita, así como a la 
observación desde cualquier punto del globo. Para ello sólo hay que modificar los 
datos de las líneas 210 a 250 y 170, respectivamente. 


MENE E MEME DIED MENTE DEI DE DEI IMEI MIME ME 
COMETA HALLEY 


1985 


kx ok o ok 
kx * xx 


VEDETTE DE DE DE METE IE DE TE DE IE DE TE DE DEDO FE DE DEDE DE METDE METE DE DEE TENE IE DE 


80 KEY OFF:*COLOR 15,1,0:WIDTH 40:CLS:DEFDBL D:DEFINT I-N 
90 D8=-1D-30: DLN=1/LOG(10+) 


100 
110 
120 
130 
140 
150 
160 
161 
162 
163 
170 
180 
190 
200 
210 
220 
224 
223 
226 
2er 
228 
227 
230 
240 
230 
260 
270 
280 
290 
300 
310 
320 
330 
340 
3530 
360 
370 
380 
390 
400 
410 
420 
4.30 


DEG=57.2957795134: DN=. 99726956642358 

DPI=3. 1415926535897934: DP12=2*DP 1 

DEF FNH(DX)=DX-24* INT (DX/24) 

DEF FND(DX)=DX-DP12*INT(DX/DPI2) 

REM 27 


RA A A o 
REM DLA=LONGI TUD 
REM DFI=LATITUD 
REM IZO=HORA DE LA ZONA 
DLA=. 2467: DF1=40. 45/DEG: 1Z0=0 


PE A A A AAA A 
DOQ=.5871023%: DEC=.9672751%+ 
DW=111.84694+*/DEG: DO=38. 1441468 /DEG 

REM DEPOCH=EPOCA 

REM DO=DISTANCIA DEL PERIHELIO 

REM DEC=ESCENTRICIDAD 

REM DW=ARGUMENTO DEL PERIHELIO 

REM DO=LONGITUD DEL NODO ASCENDENTE 

REM DI=INCLINACION 

DI=162. 239338 /DEG: DEPOCH=2446470.740627+$ 

DAC=D0Q/ (1-—DEC) 

DNC=.0172021221873%/DAC/SOR (DAC) 

REN A ARA AAA AAA TIA 


LOCATE 1,10: PRINT "COMETA HALLEY" 
LOCATE 3,3: PRINT "INTRODUCE LA FECHA Y LA HORA: " 
INPUT" ANO r ",MY 


INPUT" MES 2: ",M 

INPUT" DIA 2: ",MD 

INPUT" HORA 2 ",MH 

INPUT" MINUTO: "“,MI 

REA A 
REM PROGRAMA PRINCIPAL 

REA A A A 


DL=MH+M1/60:GOSUB 1120:PRINT 
DL=FNH (DL) : DTU=DL-1ZO 

GOSUB 1840: D4=DL : D5=0: GOSUB 870 

PRINT "FECHA :"5:PRINT USING "4H. "5MD5M5 
PRINT USING "4tHbtt. "5 MY5 | 


85 


440 PRINT TAB(20)53 "HORA 25 

450 PRINT USING "*$*"5 IH3 IM:PRINT 
460 DJ=DJ-12Z0/24+*: DL=0 

470 GOSUB 1190:D4=DA: DS5=DD: GOSUB 870 


480 PRINT "ASCENSION RECTA (h mind: "5 
4790 PRINT USING "*$%$"51H5 1M 
500 PRINT "DECLINACION to “31 


310 PRINT USING"+HH%*"5 INT C(1S)5 INT (JM) 

320 PRINT "DISTANCIA GEOCENTRICA (VA) 2 "5 
530 PRINT USING "$4+4+.$4”5DG 

540 PRINT "DISTANCIA HELIOCENTRICA (UA): ";3 
390 PRINT USING "H$HH.+HHF"5DS 


5360 PRINT "MAGNITUD TOTAL E 
970 PRINT USING "*$+H4.+"5DM1 
580 PRINT "MAGNITUD DEL NUCLEO E 
590 PRINT USING "$HH.+*"5DM2 
600 PRINT "ANG. TIERRA-SOL DEL COMETA 2: *"s 


610 PRINT USING "$%$*"35 INT (DPSI) 
620 DIJ=DJ—DTU/24+: GOSUB 1050 
630 GOSUB 1720 


640 PRINT "ALTITUD tod: ”; 
650 PRINT USING "+***"3AL 
660 PRINT "AZIMUTH (od: "5 


670 PRINT USING "*+HH*"5AZ 

680 DJ=DJ+1Z0/24+$ 

6790 IF ABS(DJ-D8)<.0001 GOTO 820 

700 GOSUB 1190 

710 PRINT:PRINT "SALIDA (h min): "; 
720 GOSUB 1940 

730 1F K=1 THEN TR=0:3 TS=24 

740 1F K=2 THEN TR=-1:2 TS=-1 

750 IF TR<O THEN PRINT " —— —-— ":GOTO 780 
760 D4=TR :D5=0:G0OSUB 870 

770 PRINT USING "$%$%*"5 IH5 INTCIM) 

780 PRINT "PUESTA (h min): “5 
790 1F TS<O THEN PRINT " -— —— ":G0TO 820 
800 D4=TS:D5=0:GOSUB 870 

810 PRINT USING "+$Hé"5IH5 INT (IM) 

3820 D8=DJ:GOTO 290 


BIO RE A 
840 REM CONVERSION DE HORAS A MINUTOS 

850 REM Y GRADOS A MINUTOS DE ARCO 

BO RE A _ _ _—_—— 


870 IH=INT (D4) : IM= (D4-1H) +60 
880 IF IM=60 THEN IM=0: IH=IH+1 


890 
900 
910 
920 
930 
940 
950 
960 


970 


980 
990 


1Z=SGN (D5) : D=ABS (DS*DEG) 

1S=INT (D) : JM= (D-15)*60 

IF JM=60 THEN JM=0: 1S=1S+1 

1S=1Z*15:2 RETURN 
O E E E 
REM ECUACION DE KEPLER 

REA A A A AA 
D=DM 

DU=D+ (DM—D+DE+*SIN(D) ) / (1—DE*COS (D) ) 


IF ABS(DU-—D)>.00001 THEN D=DU:GOTO 970 
DV=2*ATN(SOR ( (1+DE) / (1-—DE) )*TAN (DU/2)) 


1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1:270 
1.280 
1290 
1300 
1310 
1320 
1330 


RETURN 


REM HORA SIDERAL EN GREENWICH 

REM A LAS 00 hs TU 

REN AA A 
DT= (DJ-24150204$) /36523% 

D= ( . 00000107 5H$*DT+100.00213594) *DT 
D=D+.2769193989*$: 5S0=24* (D-— INT (D) ) 

RETURN 


RE, A A A A A RÁ 
IF M<3 THEN MY=MY-1:M=M+12 
DJ=1720981.5t+INT (365. Z5HxMY) 
DJ=DJ+INT (30. 60014x* (M+1)) +DL /24+MD 

RETURN 


REM 
D9=DJ-DL /244t: DT=(D9-2415020!)/36525! 
D=(.0000052796214*DT+628. 3319509044) *DT 

DLS=FND (D+4. 8816279728969744) 

D=- (. 000000057 596HHDT+. 000002517991) *DT 

DM=FND ( (D+628. 301945734) *DT+6. 25658354) 

D= (. 00000000878$+*DT-—. 00000002862344t) *DT 

DEP= (D-. 0002271114) *DT+.409319747445714 

D=-—. 0000001 264x*DT 

DE= (D+. 00004184) *DT+.016751044 

GOSUB 960 

D=FND (DLS+DV-DM> 

DRS=1 . 0000002 14tx* (1-DE*COS (DU) ) 
DXS=DRS+*COS (D) : D=DRS*SIN (D) 
DYS=D*COS (DEP) : DZS=D*SIN (DEP) 

REA O A 
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1340 REM COMETA 

LISO RE q q E 4 A 
1360 DF=COS (DO) : DG=S IN (DO)*COS (DEP) 

1370 DH=SIN(DO)*SIN(DEP) 

1380 DP=-SIN(DO)*COS(DI) 

1390 D=C0S(DO)*COS (DI) *COS (DEP) 

1400 DO=D-SIN(DI)*SIN(DEP) 

1410 D=COS (DO)*COS (D1)*SIN(DEP) 

1420 DR=D+SIN(DI)*COS(DEP) 

1430 DA2=SOR (DF*DF+DP+*DP) 

1440 DB2=SOR (DG*DG+DO*DD) 

1450 DC2=SOR (DH*DH+DR*DR) : SD=SGN (SIN (DO) ) 
1460 DA1=ATN(DF/DP) : DB1=ATN(DG/DO) 

1470 DC1=ATN(DH/DR) 

1480 IF COS(DO)<0 THEN DA1=DPI-DA1 

1490 IF SD*DB1<0 THEN DB1=DB1+DPI 

1500 IF SD*DC1<0 THEN DC1=DC1+DP1 

1510 DM=DNCx* (D9-DEPOCH) 

1520 DE=DEC:GOSUB 960 

1530 DS=DAC* (1—DEC*COS (DU) ) : DPC=DW+DV 
1540 DX=DS*DA2*SIN(DA1+DPC) 

1550 DY=DS*DB2xXSIN(DB1+DPC) 

1560 DZ=DS*DC2*SIN(DC1+DPC) 

1570 Di=DXS+DX:D2=DYS+DY: D3=DZS+DZ 

1580 DG=S0OR(D1*D1+D2*D2+D3*D3) 

1590 DA=DEG/15*ATN(D2/D1) 

1400 D=D3/DG: DD=ATN (D/SQR (1-D+*D) ) 

1610 IF D1<0 THEN DA=DA+12:GOTO 1630 
1520 IF D2<0 THEN DA=DA+24 

1430 D=.5* (DRS*DRS+DG*DG—DS*DS) /DRS/DG 
1440 DPSI=DEG*ATN(SOR (1—Dx*D) /D> 

1450 IF DPSI<O THEN DPSI=DPSI+180 

1660 DM1=5+ (S*LOG (DG)+13. 1*LOG(DS) ) *DLN 
1670 DM2=7.5+(S5*LOG (DG) +10*LOG(DS) ) *DLN 
1680 RETURN 

1690 REM === 
1700 REM COORDENADAS HORIZONTALES 
1710 REM === oo 
1720 Z=15/DEG*FNH (SO+DTU/DN-DLA-—DA) 

1730 D=COS(DFI)*C0S(DD)*C0S(Z) 

1740 D=D+SIN(DFI)*SIN(DD) 

1750 AL=DEG*ATN (D/SOR (1—D*D) ) 

1760 X=COS(Z)*SIN(DFI)-TAN(DD)*COS(DFI) 
1770 Y=SIN(Z):AZ=DEG*ATN(Y/X) 

1780 IF X<0 THEN AZ=AZ+180:GOTO 1800 


1770 1F Y<O THEN AZ=AZ+360 

1800 RETURN 

UE A A AS AN 
1820 REM DATOS DEL CALENDARIO DE DIAS JULIANOS 

LE. RE A A A A pi 
1840 D=INT (DJ+.5): A=INT ( (D/36524.3)-51. 1226) 

1850 A=D+1+A- INT (4/4) :B=4+1524 

1860 C=INT ( (B/365.25)-.3343) 2: S=INT (365. 25*C) 

1870 E=INT ( (B-S) /30. 61) 2: MD=B-S- INT (30.61x*E) 

1880 M=E-1:2:MY=C-4716:1F M>12 THEN M=M-12 

1890 IF M<3 THEN MY=MY+1 

1900 RETURN 

CTO AER Sd ARAN EA RAEE ARANA AA 
1920 REM HORARIO DE SALIDAS Y PUESTAS DEL SOL 

TO RE A AA ei 
1940 TK=-99: D7=DA: D6=DD 

1950 GOSUB 2140: IF K>0 GOTO 2090 

1960 T=DN*FNH (DA-TT-SO+DLA+1Z0) 

1970 DT=ABS(T-TK) 

1980 IF TK>0 AND DT>22 THEN T=-1:GODTO 2010 

1990 IF DT<.01*% THEN GOTO 2010 

2000 DL=T:GOSUB 1190: TK=T:GOTO 1950 

2010 TR=T: TK=-99: DA=D7 : DD=D6 

2020 GOSUB 2140: 1F K>0 THEN GOTO 1090 

2030 T=DN*XFNH (DA+TT-SO+DLA+1Z0) 

2040 DT=ABS (T—-TK) 

2050 IF TK>0 AND DT>22 THEN T=-1:GOTO 2080 

2060 IF DT<.01%* THEN GOTO 2080 

2070 DL=T:GOSUB 1190: TK=T:GOTO 2020 

2080 TS=T 

2090 RETURN 

ZOO REM AAA ii 
2110 REM HORA ANGULAR CORRESPONDIENTE AL HORARIO 
2120 REM DE SALIDA O PUESTA DEL SOL 

LLO BEN AA AA AF NAAA 
2140 CT=-.01018091142*-SIN(DD)*SIN(DFI) 

2150 ET=CT/COS (DD) /COS(DFTI):K=0 

2160 IF CT>i THEN K=2: RETURN 

2170 1F CT<-1 THEN K=1: RETURN 

2180 TT=DEG/135*ATN (SQR (1-CTCT) /CT) 

2190 IF CT<O THEN TT=TT+12 

2200 RETURN 


Apéndice C 

Adaptación de los 
programas para MSX 
a otras máquinas 


Amstrad 


Los usuarios del Amstrad podrán adaptar los programas escritos para MSX 
contenidos en este libro, siguiendo las próximas indicaciones. El principal proble- 
ma es la adaptación de la sentencia SPRITE al sistema del Amstrad, que sólo se 
acerca a este tipo de subrutinas con el uso de la orden SYMBOL. Para 
desenmarañar la lista de datos tipo CHRS() que definen a cada sprite, ha de 
emplearse la orden SYMBOL, empleando el argumento CHRS() como un dato 
más de cada símbolo. 

Para el movimiento de las figuras sobre la pantalla, el método más eficaz es 
utilizar la orden TAG. Así, cada vez que una figura tenga que desplazarse por la 
pantalla, localizaremos (LOCATE) el cursor en la nueva posición. Hay que tener 
en cuenta que cuando se mueve una figura en la pantalla no se borra la anterior. 
La mejor forma de evitar esto es definiendo cada símbolo con los pixels del borde 
actuando como ceros, según la siguiente figura. 

Cada vez que se quiera escribir texto hay que desactivar el modo de gráficos 
mediante la orden TAGOFF. 

La sentencia CIRCLE del MSX es un poco más difícil de adaptar. En este 
caso, hay que crear una subrutina que dibuje círculos excéntricos. La excentrici- 
dad se consigue variando los números que multiplican al seno y al coseno de la 
subrutina. 

Una manera eficaz de incluir a la vez texto y gráficos en pantalla sin utilizar la 


Apple ll 


orden TAGOFF consiste en definir dos ventanas; una para los textos (las cuatro 
últimas líneas) y otra para los gráficos. 
El efecto de la orden LINE del sistema MSX se puede reemplazar fácilmente 
en el Amstrad con la utilización combinada de sentencias MOVE y DRAW. 
Un programa sin gráficos como el del apéndice B para el IBM PC resultará 
muy fácil de adaptar a Amstrad. 


El Apple es uno de los primeros ordenadores personales que salieron al 
mercado y, por tanto, su sistema de gráficos no puede competir con los modernos 
ordenadores personales. El problema principal es que, al trabajar en color, la 
resolución se reduce a la mitad. Por otra parte, el Apple tampoco admite la 
utilización de sprites y hay que simularlos por medio de la creación de shapes 
tables. 

El modo de baja resolución se habilita mediante la orden GR. Este modo 
consiste en una cuadrícula de 40 x 40 pixels y deja cuatro líneas de texto que no 
borra. Tiene 16 colores que se pueden utilizar simultáneamente. 


INSTRUCCIONES EN EL MODO GRAFICO DE BAJA RESOLUCION 


COLOR = X Habilita el color X (de O a 15) para la siguiente 
instrucción. 
PLOT X, Y Sitúa un punto con el color determinado por la 


instrucción anterior en la posición de coordena- 
das X, Y. Ambas coordenadas van de O a 39, 
La posición 0,0 está situada en la parte superior 
izquierda de la pantalla. 

HLIN X1,X2 AT Y Dibuja una línea horizontal desde el punto 
X1, Y hasta el punto X2, Y. El color en el que 
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aparecerá la línea será el que se defina previa- 
mente con la orden COLOR. 


VLIN Y1, Y2 AT X Se utiliza de la misma manera que la instrucción 
anterior, pero para definir líneas verticales. 

SCRN (X, Y) Es una función que devuelve el color de punto 
referenciado. 


El modo de alta resolución se habilita mediante la orden HGR o HGR2, 
según la página que deseemos utilizar. HGR limpia la parte superior de la 
primera página de gráficos de alta resolución y la deja en negro para permitirnos 
trabajar sobre una pantalla de 280 x 160 puntos con cuatro líneas inferiores 
destinadas a texto; HGR2 limpia toda la pantalla de gráficos de alta resolución y 
la deja en negro; corresponde a la segunda página de gráficos y contiene 
280 x 192 puntos. 


INSTRUCCIONES EN EL MODO GRAFICO DE ALTA RESOLUCION 


HCOLOR = X Sitúa el color X (de O a 7) para la siguiente 
instrucción de dibujo. 
HPLOT X, Y Sitúa un punto de color previamente definido en 


el punto de coordenadas X, Y. El rango de X es 
de 0 a 279 y el de Y de 0 a 159 en HGR (de 0 a 
191 en HGR2). La posición 0,0 está situada en 
la parte superior izquierda de la pantalla. 

HPLOT X1, Y1 TO X2,Y2 Dibuja una línea desde el punto X1, Y1 hasta el 
punto X2, Y2. 

SHLOAD Carga una tabla de shapes desde casete. 

DRAW n AT X, Y Dibuja una shape número n en la posición X, Y 
en el color determinado en la última orden 
HCOLOR. 

* XDRAW n AT X, Y Dibuja una shape, la n, en el punto X, Y con el 
color complementario al determinado en la or- 
den HCOLOR. 

ROT = X Sitúa la rotación de la shape para las órdenes 
DRAW y XDRAW. El rango de X es de O a 
255. X=0 es vertical; los demás dependen de la 
escala. 

SCALE = X Sitúa la escala de 1 a 255 para el DRAW y el 
XDRAW. 


Una shape es una figura que se dibuja en la pantalla de alta resolución 
guiándose por una tabla de bytes en memoria. En este apéndice no nos detendre- 
mos en la creación de estas figuras, por lo cual remitimos a los usuarios de Apple 
a la lectura del capítulo 9 del manual de referencia de programación en BASIC 
Applesoft que acompaña al ordenador. 


Spectrum 


La adaptación del BASIC MSX al BASIC de Sinclair en los programas de este 
libro es, básicamente, un trabajo de modificación del sistema gráfico de represen- 
tación de sprites. La potencia del estándar MSX para representar figuras en 
movimiento en pantalla no es fácil de implementar en una máquina como el 
Spectrum. Sin embargo, los dibujos de los programas del Halley son lo suficiente- 
mente sencillos como para poder ser representados con los parcos comandos del 
BASIC Sinclair. Los dibujos en alta resolución estáticos no presentarán problema 
al programador medianamente experimentado. La mejor referencia para adaptar 
los programas escritos para el sistema MSX al BASIC del Spectrum es estudiar 
las diferencias entre las dos versiones escritas para estos ordenadores de los 
programas “Halley” y “El Halley sobre el cielo nocturn1o” (4 y 5, respectivamente). 

Los comandos PRESET y PSET de MSX trazan un punto en coordenadas 
absolutas en color de papel y tinta, respectivamente; por tanto, son fácilmente 
sustituibles mediante el PLOT del Spectrum. LINE es equivalente a DRAW en 
BASIC Sinclair, pero puede aparecer con coordenadas relativas o absolutas, según 
haya uno o dos pares de parámetros. Si aparece con una pareja de parámetros se 
puede sustituir directamente LINE por DRAW, teniendo en cuenta únicamente 
que en MSX se dispone de 191 puntos verticales. Si aparecen dos pares de 
parámetros, el primero indica la posición inicial y el segundo la final. El comando 
DRAW de MSX se emplea con el macrolenguaje gráfico, no existente en el 
Spectrum. Por su complicación se recomienda la consulta de algún manual para 
su comprensión (véase Guía del programador y manual de referencia MSX, de esta 
misma editorial). Por último, CIRCLE es capaz de trazar círculos y elipses en MSX, 
por lo que, en este último caso, deberá generarse la elipse en el Spectrum punto a 
punto mediante su ecuación. 

Los sprites del sistema MSX pueden sustituirse mediante los gráficos definidos 
por el usuario del Spectrum. El movimiento de éstos debe implementarse median- 
te impresión y borrado alternativos, procurando que el tiempo que transcurra 
entre un borrado y una nueva reimpresión sea mínimo. Si no se usa una rutina 
especial de las muchas aparecidas en distintas publicaciones, la impresión de las 
figuras deberá limitarse a una de las posibles posiciones de baja resolución del 
Spectrum. La calidad obtenida no será la misma, pero la solución es válida. 

El resto de las adaptaciones necesarias no debe presentar gran dificultad. Aquí 
podemos nombrar algunas de las más notorias como, por ejemplo, la necesidad 
de anteceder las sentencias de asignación con un LET. La sentencia DIM para 
variables alfanuméricas varía en el Spectrum en lo relativo al último parámetro, 
que indica la longitud de cada cadena. Haciendo ésta suficientemente grande no se 
planteará ningún problema. El USING de MSX debe implementarse en Spectrum 
mediante TAB y AT. ON GOTO y ON GOSUB no existe tampoco en el BASIC 
Sinclair, pero pueden sustituirse fácilmente con IF. Podríamos seguir nombrando 
gran cantidad de pequeños detalles, pero la parquedad de este pequeño libro nos 
lo impide. Para cualquier duda que surja en la adaptación es recomendable 
recurrir a un buen manual de MSX, como el nombrado anteriormente. 
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Apéndice D 
Programas comerciales 
de astronomía 


Hay programas comerciales de astronomía diseñados para la comprensión de 
los conceptos astronómicos. Son los programas educacionales (E). También hay 
otro tipo de programas de tipo “planetarium” (P) que nos permiten observar toda 
la bóveda celeste en la pantalla del ordenador, localizando con sus coordenadas 
exactas los planetas, cometas y constelaciones en una fecha indicada. Hasta te 
mostrarán la Luna en el tamaño exacto para la fecha, día y hora indicada y su 
altura sobre o bajo el horizonte. Con ello obtendrás un mapa celeste completo a 
la vez que un calendario astronómico con miles de años de duración. Por último 
hay programas tipo, simulación o videojuego (J). 


PROGRAMA/ORDENADOR TIPO PRODUCIDO POR 
STARWARE E Atari Program Exchange 
Atari 400-800 Box 3705 
Santa Clara, CA 95055 
EE.UU 
SELLAR ASTRONOMY E Cross Educational Software 
/Apple II+, IBM Box 1536 


Ruston, LA 71270 
EE.UU 


SOLAR SYSTEM ASTRONOMY 
Apple Il, IBM PC 


THE OBSERVATORY 
[Apple 1I+, Ile 


SKY TRAVEL 
[Commodore 64 


TELL STAR 
Apple II, IBM PC, 
Pejr, Macintosh 


THE HALLEY PROYECT 
Commodore 64 
Attari 


Cross Educational Software 
Box 1536 
Ruston, LA 71270 
EE.UU 


Lightspeed Software 
2124 Kittredge St Suite 185 
Berkeley, CA 94704 
EE.UU 


Commodore 
1200 Wilson Dr. 
Westchester, PA 19380 
EE.UU 


HoloByte Corp. 
1050 Walnut, Suite 325 
Boulder, CO 80302 
EE.UU 


Mindscape 
3444 Dundee Rd. 
Northbrook, IL 60062 
EE.UU 
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El mayor acontecimiento de los ochenta es la nueva visita 
del cometa HALLEY. 


En esta ocasión el cometa se encuentra con una generación capaz 
de analizarle de forma sofisticada mediante el uso de ordenadores 
personales. 


COMETAS EN TU MICRO te da todos los programas necesarios 
para que puedas predecir y analizar el comportamiento 
del HALLEY o de cualquier otro cometa. 


La primera parte del libro ofrece una panorámica general sobre 
la estructura del sistema solar, la naturaleza de los cometas 
y datos sobre el HALLEY y su observación práctica. 


En la segunda parte encontrarás cinco programas gráficos | de. 
y dos numéricos con los que analizar exhaustivamente el paso A 
del HALLEY o de otros cometas. " 


. Si te atrae la astronomía cometaria y posees un microordenador, 
, con COMETAS EN TU MICRO profundizarás en el diseño 
de programas astronómicos en BASIC y podrás comprender mejor 
las leyes que rigen el sistema solar y los cometas. 
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